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ABSTRACT 

Cervantes  de  Gortari,  Jaime  Gonzalo.  Ph.D.,  Purdue 
University,  August  19  76.  A'J  EXPERIMENTAL  STUDY  OF  THE 
FLAPPING  MOTION  OF  A TURBULENT  PLANE  JET.  Major  Professor: 
Victor  W.  Goldschmidt. 

Measurements  were  conducted  in  order  to  characterize 
the  flapping  motion  of  a turbulent  plane  jet.  The  flapping 
is  defined  as  the  (apparently  natural)  lateral  pseudo- 
oscillatory  motion  (in  the  average)  of  the  flow  field  of  • 
the  jet.  The  technique  employed  in  these  measurements  was 
to  compute  crosscorrelation  functions  between  the  velocity 
fluctuations  (one  delayed  in  time)  in  the  longitudinal  di- 
rection, at  two  points  on  opposite  sides  of  the  jet.  Stand- 
ard hot-wire  anemometry  and  on-line  digital  processing  in- 
strumentation were  used.  Four  types  of  measurements  v/ere 
performed:  (a)  with  the  probes  symmetrically  positioned 

with  respect  to  the  centerline  of  the  jet;  (b)  with  the 
probes  asymmetrically  set  with  respect  to  the  centerline; 

(c)  with  the  probes  at  points  of  different  longitudinal  co- 
ordinate; (d)  with  an  x-wire  at  the  centerline  of  the  jet 
where  the  autocorrelation  of  the  lateral  component  of  the 
velocity  was  measured. 

The  existence  of  flapping  was  indicated  by  a dis- 
tincti\  negative  value  at  zero  time  delay  and  a quasi- 
periodic  nature  at  the  measured  crosscorrelation  functions. 


XV 


These  were  obtained  at  longitudinal  stations  ranging  from 
10  to  100  times  the  slot  width.  The  frequency  characteriz- 
ing the  flapping  motion  showed  approximate  self-preservation 
(when  scaled  by  the  half-width  and  the  mean  velocity  at  the 
centerline  of  the  jot)  for  stations  with  longitudinal  co- 
ordinate larger  than  30  times  the  slot  width.  It  was  found 
that  fb/U^  0.11.  The  flapping  frequency  did  not  show  de- 
pendence on  the  lateral  coordinate.  The  amplitude  of  flap- 
ping was  estimated  to  about  20%  of  the  jot  half-width.  The 
flapping  behavior  seems  to  travel  in  the  downstream  direction 
with  a velocity  smaller  than  the  convective  velocity  of  the 
turbulent  structure.  The  frequency  of  flapping  showed 
Reynolds  number  independence  in  the  range  7900  £ Re  £ 15100, 
where  Re  is  based  on  the  nozzle  exit  characteristics. 
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CHAPTER  1 
INTRODUCTION 

The  subject  of  turbulent  shear  flows  has  received 
considerable  attention  in  recent  years.  Sufficient 
motivation  for  studying  this  kind  of  flow  will  result  if 
one  recalls  engineering  applications  like  mixing  processes, 
diffusers,  jet  engines,  fluid  amplifiers  or  the  more  recent 
attempts  to  control  environmental  pollution,  to  mention  but 
a few. 

The  understanding  of  turbulence  in  general  and 
turbulent  shear  flows  in  particular,  is,  however,  far  from 
complete.  The  main  reason  for  this  is  the  lack  of  a general 
mathematical  theory  of  turbulence.  The  advances  in  this 
field  have  proceeded  to  a large  extent  on  an  experimental 
basis.  By  this  means,  empirical  models  are  constructed 
which  may  be  used  in  a descriptive  way  or  as  a criterion  for 
design.  Experimental  results  also  contribute  to  the  general 
body  of  data  against  which  any  mathematical  theory  may  be 
verified. 

Among  the  various  aspects  of  turbulent  shear  flows 


not  sufficiently  investigated  up  to  date  is  the  phenomenon 
of  flapping  of  a turbulent  plane  jet.  The  flapping  is 
defined  as  the  lateral  oscillatory-like  motion  (in  the 
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averaye)  of  the  flow  field  of  the  jet.  Kffectivcly,  a plane 
jet  issueing  from  a rectangular  slot  attains  a turbulent 
state  after  some  distance  downstream  from  the  nozzle  exit. 

This  turbulent  flow  is  characterized  by  random  fluctuations 
in  the  velocity  field  and  the  usual  statistical  parameters 
like  moan  velocity,  intensity  of  turbulence,  energy  spectra 
and  length  scales  can  be  defined  and  measured.  However, 
hidden  in  the  randomness  of  the  longitudinal  velocity 
fluctuations  appears  to  exist  a low-f roquency , almost 
periodic  component.  The  longitudinal  velocities  detected  at 
two  points,  each  one  at  opposite  sides  of  the  jet  centerline, 
are  such  that  their  corresponding  low-frequency  components 
are  in  counterphase  with  respect  to  each  other  (in  the 
average) . In  other  words,  when  the  velocity  at  one  point 
increases,  the  velocity  at  the  opposite  point  generally 
decreases.  The  overall  result  is  that  the  flow  field  moves 
in  a lateral  oscillatory  fashion,  as  depicted  in  Figure  1.1. 

The  objective  of  this  investigation  was  to  perform  experi- 
ments aiming  to  detect  and  quantify  the  phenomenon  of  flapping 
of  a turbulent  plane  jet.  The  technique  employed  in  this 
investigation  was  to  compute  statistical  crosscorrelations 
between  longitudinal  components  of  the  velocities  (one 
delayed  in  time)  at  two  points,  each  one  on  opposite  sides 
of  the  jet.  The  existence  of  flapping  would  be  indicated  by 
a distinctive  negative  valvic  at  zero  time  delay  and  a 
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iping  of  a Turbulent  Plane 
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quaai-periodic  risiturc  of  the  moasurocl  croascorrolation 
functions.  Standard  hot-wire  anemometry  and  on  lino  digital 
processing  instrvimontation  were  used  in  the  four  typos  of 
performed  measurements:  (a)  with  two  normal  hot-wires 
symmetrically  positioned  with  respect  to  the  centerline  of 
the  jet;  (b)  with  the  probes  asymmetrically  sot  with  respect 
to  the  centerline;  (c)  with  tlie  probes  at  points  of 
different  longitudinal  coordinate;  (d)  with  an  x-wire  at  the 
centerline  of  the  jet  where  the  autocorrelation  of  the 
lateral  component  of  the  velocity  was  measured.  Longitudinal 
stations  ranging  from  10  to  100  times  the  slot  width  were 
investigated. 

These  experiments  were  planned  to  satisfy  the 
following  specific  objectives.  First,  to  clearly  establish 
the  existence  of  the  flapping  motion  as  opposed  to  a puffing 
effect,  i.e.,  longitudinal  oscillations  of  the  flow  field 
(this  effect  would  be  originated  by  artificial  instabilities 
somewhere  upstream  of  the  nozzle  exit  and  would  bo 
characterized  by  periodic  crosscorrelation  functions  with 
positive  values  at  zero  time  delay) . Second,  to  look  for 
any  dependence  of  tlio  flapping  frequency  on  the  longitudinal 
and  lateral  coordinates.  Third,  to  study  the  traveling 
characteristics  of  the  flapping  behavior  in  the  downstream 
direction.  Finally,  to  show  tlie  possibility  of  Reynolds 
number  independence  and  self-preservation  of  tlie  flapping 
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motion . 


In  what  follows,  after  a review  of  the  general 
characteristics  of  turbulent  plane  jots  and  relevant 
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correlation  techniques,  a summary  of  the  available 
literature  related  to  this  investigation  is  presented. 
Thereafter  a detailed  description  is  given  of  the  experimental 
procedures  and  employed  instrumentation.  The  last  three 
chapters  refer  to  the  experimentally  obtained  results  and 
some  concluding  remarks  are  offered.  Suggestions  for 
further  work,  in  the  light  of  this  investigation,  are  also 
given . 

It  is  hoped  that  characterizing  the  flapping  behavior 
of  turbulent  plane  jets  may  contribute  to  a better  insight 
of  the  structural  properties  of  turbulence,  particularly  in 
the  context  of  the  new  trends  in  turbulence  research,  where 
a more  deterministic  point  of  view  is  employed. 


CHAPTER  2 

REVIEW  OF  CHARACTERISTICS  OF  PLANE  JETS 
AND  CORRELATION  TECHNIQUES 

2 . 1 General  Characteristics  of  a Plano  Turbulent  Jet. 

A plane  turbulent  jet  belongs  to  the  class  of  flows 
known  as  free  turbulent  shear  flows.  Typically,  a plane 
jet  is  formed  when  a fluid  with  sufficiently  large  velocity, 
is  discharged  through  a slot-shaped  orifice  into  a volume  of 
quiescent  fluid.  Initially,  a wedge-shaped  core  is  produced 
close  to  the  nozzle  with  a uniform  velocity  and  low  turbu- 
lence level,  bounded  by  two  mixing  layers  with  high  turbu- 
lence intensity,  which  gradually  diffuse  momentum  from  the 
core  to  the  surroundings  (see  Figure  2.1). 

As  the  development  continues,  the  two  mixing  layers 
merge  together  after  a distance  of  about  six  times  the  slot 
width  and  the  mean  velocity  at  the  centerline  starts  to  de- 
cay. After  this  transitional  region  a fully  turbulent  re- 
gion is  established,  where  the  mean  velocity  profiles  take  a 
shape  close  to  a gaussian  or  normal  function  and  the  static 
pressure  is  constant. 

In  the  fully  turbulent  region,  the  jet  attains 
Reynolds  number  similarity,  i.e.,  the  mean  flow  structure 
is  independent  of  viscosity  and  the  transverse  distributions 
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of  mean  velocity  and  other  mean  quantities  (although  chang- 
ing in  the  downstream  direction)  become  self-preserving  if 
they  are  expressed  non-dimensionally  through  prop  scales. 
For  instance,  the  mean  velocity  profile  becomes  an  universal 
function  if  (a)  the  mean  velocity  itself  is  divided  by  the 
jet  centerline  mean  velocity  and  (b)  the  corresponding 
transverse  y coordinate  by  the  jet  half-width  b.  The 
latter  is  defined  as  the  lateral  distance  where  the  mean 
velocity  is  exactly  one-half  of  the  centerline  value  (see 
Figure  2.1). 

The  basic  equations  of  motion  can  be  simplified  for 
this  kind  of  flow  by  making  use  of  one  experimental  observa- 
tion: the  flow  is  confined  to  a region  of  space  elongated 

in  the  direction  of  the  general  motion  and  narrow  in  the 
transverse  direction.  In  effect,  through  an  order  of  magni- 
tude analysis  (Townsend  [1956]),  the  incompressible  mean 
flow  equations  can  be  reduced  to: 


momentum  in 
the  mean  flow 
direction: 


U 

9x 


9 (uv) 
9y 


0 


(2.1) 


continuity:  ^ = 0 (2.2) 

dx  dy 

Similarity  solutions  of  equations  (2.1)  and  (2.2) 
were  first  obtained  by  Tollmien  [1926],  as  referred  by 
Schlichting  [1968],  after  using  Prandtl's  mixing  length 
hypothesis.  Gocrtler  [1942]  and  Reichardt  [1951]  also 
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where  n = 7-67  b is  the  half-width  of  the  jet  (Fiejuro 

2.1)  and  Uj,^  is  the  mean  velocity  at  the  centerline.  Equa- 
tion (2.4)  is  plotted  in  Appendix  A toqether  with  measure- 
ments of  U taken  in  the  present  jet. 

Further  details  about  self-preservation,  basic  equa- 
tions and  various  other  aspects  of  turbulent  jets  are  avail- 
able in  the  literature  (Townsend  (1956),  Abramovich  [19631, 
Schlichting  (19681). 


2 2 Correlation  Techniques 

Correlation  analysis  is  a well  accepted  statistical 
technique.  Its  usage  in  sciences  like  agriculture,  biology 
and  economic  statistics,  (altliough  somewhat  restricted 
through  the  concept  of  correlation  coefficient),  is  exten- 
j sively  used  in  comparing  series  of  measurements  of  any  kind. 

li 

It  is  perliaps  in  the  theory  and  applications  of  communica- 
tions and  random  processes  where  its  broadest  utilization 
may  bo  found.  (see  for  instance  Lange  [19661,  Lathi  [1968], 
|(  Bondat  and  Piorsol  [19711,  Jenkins  and  IVatts  [19681). 

j Correlation  techniciues  liave  been  also  widely  used  in 


J 


o.xper imontal  turbulent  flow  studies  since  Taylor  [1935) 
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first  suggested  them.  Indeed,  one  of  the  methods  used  to 
get  some  insight  into  the  structural  and  transport  char- 
acteristics of  this  kind  of  flows  is  to  measure  correlations 
between  various  fluctuating  quantities  like  velocity, 
pressure,  temperature,  etc.  at  different  points  and  times. 

The  correlation  function  gives  the  general  relation- 
ship between  two  random  variables.  In  a mathematical  sense 
it  is  defined  as  follows. 

Consider  a random  process  {x(t)}  as  described  by  its 
ensemble  and  a probability  measure  over  the  ensemble.  Fur- 
ther, consider  two  random  variables  and  X2  contained  in 
the  process  and  their  joint  probability-density  function 
p(Xj^,X2).  These  random  variables  x^^  = x(t^)  and  X2  = ^^^2^ 
are  defined  by  the  amplitudes  of  the  sample  functions  form- 
ing the  ensemble,  at  instants  t^^  and  t2,  respectively.  The 
correlation  between  the  two  random  variables  x^  and  X2,  or 
better,  the  autocorrelation  of  the  random  process  {x(t)}  is 
defined  as  the  ensemble  average: 


Rx(T) 


.00 

— 00 


CO 

Xj^X2P  (Xj^,X2)  dx^dx2 

— 00 


(2.5) 


where  t is  the  difference  between  t^  and  t2. 

If  the  random  process  is  ergodic,  the  complete  sta- 
tistical properties  of  the  process  may  be  described  through 
the  properties  of  any  one  of  the  sample  functions.  In  other 
words,  averaging  over  the  whole  ensemble  is  equal  to 
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timo-avoraqing  over  a single  sample  function.  The  autocorre- 
lation function  then  becomes: 


f 


% 


= lim  if 


X (t) X (t+x ) dt 


-T 


Vo 


(2.6) 


where  x(t)  is  the  amplitude  at  any  instant  t for  a single 
sample  function. 

The  above  concepts  may  be  extended  to  the  case  of 
two  jointly  ergodic  random  processes  {x(t)}  and  {y(t)}.  The 
crosscorrelation  between  the  two  processes  may  be  defined 
in  terms  of  the  time-average: 


1^2 


R ( T ) = 1 im 
xy 


x(t)  y(t+T)dt 


(2.7) 


-T 


4 


where  now  x(t)  and  y(t)  represent  any  samples  of  the 
processes  {x(t)}  and  {y(t)},  respectively. 

The  main  properties  of  correlation  functions  are 
summarized  in  the  following  relationships: 


R^(t)  = R^(-x) 


Rj^(O)  = x' 


= R^(T+nT),  for  x(t)  = x(t+nT) 


(2.8) 

(2.9) 

(2.10) 


(2.11) 
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^2.12) 

|Kjjy(T)  I 1 [Rj^(0)Ry(0)  l’’  (2.13) 

In  utilizing  the  defined  correlation  functions  in 
the  field  of  turbulence,  it  is  convenient  to  consider  the 
sample  functions  x(t)  and  y(t)  as  any  two  quantities  of  in- 
terest. For  instance,  assume  that  the  velocity  fluctuations 
at  two  points  separated  by  a distance  r are  considered; 
u(x)  and  u{x+r) , respectively,  where  x is  a space  coordi- 
i\ate;  then  the  autocorrelation  of  u(x)  is  written  as: 

R(t)  = u(x,  t) u(x, t+xT  (2.14) 

where  the  overbar  denotes  time-averaging.  The  space- 
correlation  between  u(x)  and  u(x+r)  is  defined  as: 

R(r)  = u (x, t) u (r+r , t)  (2.15) 

Furthermore,  if  both  space  and  time  are  taken  into 
account,  the  space-time  correlation  function  may  be  defined: 

R(r,T)  = u (x , t) u (x+r , t+T ) (2.16) 

Note  that  the  above  defined  space- time  correlation 
becomes  an  autocorrelation  when  the  separation  distance  is 
zero,  and  becomes  a space-correlation  for  null  time  delay  t 
(stationarity  and  homogeneity  are  assumed) , 

Normalized  (dimensionless)  correlation  functions  are 
obtained  after  dividing  them  by  the  root-moan-S(]uare  values 
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of  the  corresponding  fluctuations  (in  this  work  a normalized 
correlation  function  will  be  indicated  by  an  underlined  R) . 

A review  on  correlation  functions,  including  measurements  in 
various  flow  situations,  is  given  by  Favre  [1965]. 

Correlation  techniques  have  become  a powerful  tool 
in  studying  turbulent  flow  properties.  They  can  give  an  in- 
dication of  the  changing  pattern,  both  in  space  and  time,  of 
a given  flow.  They  can  provide  a means  of  testing  theoreti- 
cal predictions,  as  in  some  instances  of  isotropic  turbu- 
lence (Hinze  [1959]).  The  use  of  correlation  functions  in 
characterizing  the  size  and  shape  of  turbulent  eddies,  and 
estimating  appropriate  scales  for  the  corresponding  flow,  is 
generalized.  The  study  of  convective  velocities  (the  vel- 
ocity at  which  the  turbulent  structure  moves)  and  the  valid- 
ity of  Taylor's  hypothesis  (that  turbulence  remains  as  a 
frozen  strvicture  moving  at  U)  has  been  carried  out  through 
the  measurement  of  correlation  functions  for  different  flow 
configurations  (Favre  ct  al.  [1958],  Fisher  and  Davies 
[1963])  as  noted  by  Ott  [1972]. 

Finally,  the  mean  properties  of  the  large-scale 
orderly  structures  in  the  mixing  region  of  round  jets,  have 
been  investigated  using  correlation  techniques  (Lau  et  al. 
[1972],  Fuchs  [1972],  Ko  and  Davies  [1975]). 

The  use  of  correlation  tochni(]ucs  in  this  investiga- 
tion is  circumscribed  by  the  following  normalized  cross- 
correlation  function; 
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R(T)  = + 

/ u^(y)  / u^(-v) 


(2.17) 


that  is,  a time  correlation  function  between -‘xial  velocity 
fluctuations  at  opposite  sides  of  the  jet.  A positive  value 
of  R for  zero  time  delay  between  signals  (t-0) , would  indi- 
cate that  both  velocity  fluctuations  change  in  the  same 
direction  (in  the  average) , whereas  a negative  value  would 
indicate  that  when  one  fluctuation  increases,  the  other  one 
generally  decreases. 

The  same  correlation  function  was  employed  by 
Wygnans)<i  and  Gutmark  [1971]  in  trying  to  relate  the  motion 
of  the  boundaries  of  a plane  jet.  They  used,  however,  inter- 
mittency  signals  consisting  of  constant  positive  or  negative 
voltages  representing  turbulent  or  potential  flow,  respec- 
tively, in  place  of  the  velocity  fluctuations  in  equation 
(2.17).  Their  conclusion  (they  found  a zero  correlation) 
was  that  the  boundaries  of  the  jet  flow  i:iov<'d  independently 
of  each  other. 

On  the  other  hand,  Goldschmidt  and  Hradshaw  (1973) 
performed  preliminary  experiments  following  equation  (2.17) 
exactly  and  found  that  for  large  averaging  limer.,  negative 
correlation  values  at  zero  time  delay  are  per ;;  i t ent . 

The  experiments  reported  in  this  woi  k ui  c'  a natural 
extension  of  the  latter  preliminary  inve;,t  i .|,it  ion . 
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CHAPTER  3 ‘ I 

LITERATURE  REVIEW  j 


In  this  chaptor  a review  of  the  literature  most 
closely  related  to  this  investigation  is  presented.  The 
flapping  motion  of  a turbulent  plane  jet,  a naturally 
occurring  phenomenon,  has  not  been  previously  reported. 
Indeed,  a carefully  conducted  search  in  the  literature  did 
not  show,  except  for  the  preliminary  experiments  of 
Goldschmidt  and  Bradshaw  [1973] , any  reference  to  this 
aspect  of  plane  jets.  It  seems  therefore,  that  wo  are  faced 
with  an  unknown  effect  not  sufficiently  investigated  vip  to 
now. 

With  this  in  mind,  a review  of  the  literature  was 
done  along  three  directions.  One  of  them,  as  summarised  in 
the  first  section  of  this  chapter,  really  refers  to  a 
documentation  of  turbulent  plane  jots.  As  will  bo  noted, 
experimental  investigation  on  plane  jets  has  been  conducted 
to  a largo  extent  and  the  measurement  of  several  statistical 
quantities  has  boon  successfully  accomplished. 

The  other  two  areas  of  the  explored  literature  are 
the  recent  work  on  coherent  structures  and  the  study  of 
disturbed  plane  jots.  They  arc  relatively  now  linos  of 
approach  in  studying  turbulent  .shear  flows  <ind  are  presented 
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^ hero  for  tho  purpose  of  providiny  a backyround  to  the 

problem. 

In  the  appendix  a review  of  the  wave-yuido  modeliny 
of  turbulent  shear  flows  is  also  presented. 

3 . 1 Plane  Turbtilent  Jots 

The  documentation  of  plane  turbulent  jets  has  been 
cTuite  extensive.  In  what  follows,  only  ‘jomo  of  the  salient 
publications  arc  retraced. 

Classical  publication;-,  include  Townsend  [l9G5]  and 
Abramovich  [1963].  They  emphasize  the  statistical  descrip- 
tion of  the  flow  variables,  tlie  similarity  and  self- 
preservation  characteristics  and  the  role  of  phei^omenoloy i- 
I cal  theories  for  various  turbvilent  shear  flows. 

Van  der  lleynn  Zi jnen  [1958],  after  having  reviewed 
earlier  theories  and  experiments,  estimatetl  turbulent  shear 
stresses  and  the  momentum  exchange  coefficient  and  mixing 
length  in  a plane  jot. 

More  I'cccnt  experimental  results  were  reported  by 
Heskostad  [l965]  , Bradbury  [l965]  and  Mih  and  lioopes  [l972]. 
The  measurements  by  Heskostad  included  distributions  of  mean 
velocity,  turbulent  intensities,  Reynolds  stresses  and  inter 
mittoncy  and  flatness  factors  for  tlie  streamwise  velocity 
component.  An  energy  balance  for  a location  in  the  self- 
preserving  region,  was  also  done.  Bradbury  did  similar 
measuremonts  including  some  spectral  data  vm.d  compared  them 
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with  proviovuiLy  repotted  results,  particularly  the  ones  by 
Townsend  [l056]  for  a i)lane  wake.  lie  also  ronorted.  on  the 
eddy  structvire  and  the  applicability  of  .iv.iiKiblo  phenomeno- 
logical theories.  Mih  and  Hoopcs  perfc)rmeil  similar  experi- 
ments in  a water  jet.  They’  useil  an  elt'ct. I'omaijnetic  induc- 
tion method  to  detect  the  fluctuatiruj  velocities.  Their 
measurements  included  enenjy  spectra  and  length  mic:ro-  and 
macro-scales . 

Enervjy  spectra  measurements  in  a plane  jot  have  boon 
also  reported  by  Bashir  and  Uberoi  [rj7l]  and  by  Young 
[^1973j.  The  latter  also  determined  micro-  and  macro-scale 
distributions. 

The  measurement  of  static  I'rcssuri'  fluctuations  has 
been  of  interest  for  a long  time  in  turbulence  studies.  Re- 
cent attempts  to  deal  with  this  aspect  of  shear  flows  are 
reported  by  Siddon  [1069],  Spencer  [l970j  and  Planclwn 
[1974]  . 

Very  recently,  Gutmark  and  Wynnanski  [l976l  published 
their  measurements  in  a plane  jot.  They  included  mean  vel- 
ocity profiles,  turbvilent  intensities,  third  atid  fourth  or- 
der moment  quantities,  intermittency  factor  and  two- 
point  correlations.  Extensive  vise  of  conditional  sampling 
techniques  was  done.  The  resulted  data  obtained  exclusively 
within  the  turbulent  zone  of  tlie  jot  wore  compared  with  con- 
ventionally averaged  measurements. 

Table  3.1  summarizes  the  measurements  reviewed  in 
this  section.  As  noted,  the  quantities  that  have  been 


Table  3.1  Measurements  performed  in  Turbulent  Plane  Jets. 
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Other  works  not  specifically  referred  to  here  are  listed  in  the  General 
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measured  in  a turbulent  plane  jet  are:  mean  velocity  pro- 
files, turbulent  intensities,  turbulent  shear  stresses,  in- 
termittency  factors,  energy  spectra,  length  micro-  and 
macro-scales  and  some  higli-order  statistical  moments.  A 
possible  closure  theory  of  turbulence  requires  nevertheless, 
more  experimental  information;  particularly  in  respect  to 
the  structviral  character  ist  ics  of  the  turbulence  in  shear 
flows.  The  stvidy  of  the  flapping  of  a plane  jet  might  con- 
tribute to  this. 

3 . 2 Orgg ni::ed  .Uructures  in  Free  Tvirbulent  Shear  Flows 

A new  line  of  work  in  understanding  turbulent  shear 
flows  has  emerged  in  the  last  ten  or  so  years.  It  consists 
of  studying  the  large  scale  coherent  structures  that  seem  to 
exist.  The  dove lopr.ient  of  digital  data  processing  tech- 
niques and  the  utilization  of  more  sophisticated  experimen- 
tal methods,  among  other  things,  have  contributed  to  the 
growth  of  this  new  trend  in  turbulence  research.  The  in- 
terest in  these  large  scale  structures  has  been  highly  mo- 
tivated by  the  aerodynamic  noise  problem.  Indeed,  a large 
amount  of  the  available  literature  refers  to  research  con- 
ducted in  the  near-field  of  a circular  jet  and  its  relation 
to  noise  production.  Since  the  pioneering  works  by  Townsend 
[1956],  Grant  [1953]  and  Corrsin  and  Kistler  [1054],  there 
is  an  increasing  evidence  that  the  large  scale  strvictviros 
greatly  contribute  to  the  transport  of  momentum,  mass  and 
heat  in  shear  flows. 
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Among  the  flow  visualization  experiments  which  show 
how  these  ordered  structures  might  look,  are  the  now  classic 
works  by  Brown  and  Roshko  (1971]  and  Roshko  (1975).  They 
took  shadowgraph  pictures  of  a turbulent  mixing  layer  be- 
tween two  different  gases  moving  at  high  speed.  Coherent 
vortex-like  structures  appeared  in  their  pictures,  with  size 
and  space  between  them  increasing  with  distance  downstream 
of  the  origin.  The  merging  of  neighboring  eddies  accounted 
for  the  increase  of  eddy  size  and  thus  the  spreading  of  the 
mixing  layer,  although  this  coalescing  process  occurred 
quite  irregularly. 

Winant  and  Browand  [1974]  reported  similar  experi- 
ments conducted  in  a channel  where  two  streams  of  water  at 
different  velocity  were  brought  together.  Although  the 
Reynolds  number  range  was  considerably  lower  than  the  one 
used  by  Brown  and  Roshko,  again  the  coalescing  mechanism 
shown  by  dye  injection  seemed  to  control  the  growth  of  the 
mixing  layer.  They  studied  the  properties  of  these  struc- 
tures by  performing  ensemble  averages  through  conditional 
sampling,  over  a largo  number  of  structures.  They  found  a 
characteristic  Strouhal  number  fh/U  of  0.2,  associated  with 
the  passage  of  the  structures  (f  is  the  frequency  of  the 
lowest  peak  in  the  longitudinal  velocity  spectrum,  h is  the 
mixing  layer  thickness  and  U is  the  moan  value  of  tlie  two 
streams  velocities) . They  further  proposed  a model  for  the 
shear  layer  growth  based  on  the  merging  process. 


In  conducting  a statistical  study  of  tlio  pressure 
and  velocity  fluctuations  in  a mixing  layer,  dpeneer  and 
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Jones  [1971]  using  modern  signal  processing  techniques, 
noted  the  existence  of  peaks  in  the  one-dimensional  energy 
spectra  of  velocity  fluctuations.  They  found  tliese  at  the 
irrotational  outer  regions  of  the  shear  layer.  As  sugges- 
ted by  Winant  and  Browand  [1974],  these  peaks  may  be  related 
to  the  passage  of  the  paired  structures,  with  a Strouhal 
number  f U of  0.19. 

Yule  et  al.  [1974]  and  Bruun  and  Yule  [1974],  as 
described  by  Davies  [1975],  investigated  tlie  coherent  struc- 
tures in  the  first  ten  diameters  of  round  jets.  Visualiza- 
tion techniques  showed  interaction  and  coalescence  of  large 
ring  vorte.\-like  eddies  with  stronger  three  dimensional 
characteristics  than  the  previously  reported  mixing  layers. 
Measurements  of  turbulent  energy,  sliear  stresses  and  cross- 
correlations associated  with  those  structures,  wore  done. 

The  phenomenon  of  vortex  coalescence  in  various 
shear  flows  have  been  reported  previously  in  the  literature. 
Outstanding  publications  are  the  ones  by  Rockwell  and 
Niccolls  [1972]  for  a planar  jet  and  Becker  and  Mass aro  [1968] 
for  a circular  jet.  Rockwell  and  Niccolls  reported  visuali- 
zation experiments  using  the  hydrogen  bubble  technique,  in 
studying  the  breakdown  process  of  a plane  jet  witli  exit 
Reynolds  numbers  ranging  from  2,000  to  10,000.  Two  types 


of  vortex  growth  and  coalescence  were  found:  symmetrical 
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and  asymme tr ica 1 with  respect  to  the  jot  centerline.  The 
Strouhal  number  corresponding  to  the  formation  of  the  vor- 
tices was  estimated  to  follow  the  relationship  St  = 0.012 v^Rc  , 
where  both  the  Strouhal  number  St  and  the  Reynolds  number  Re 
are  based  on  the  nozzle  exit  characteristics . 

Becker  and  Massaro  [1968]  performed  experiments  with 
a circular  jot  at  moderate  Reynolds  nun±)ers.  Using  light 
scatter  and  photography  by  sheet  illumination  techniques, 
they  were  able  to  study  the  vortex  evolution  and  coalescence 
in  the  instability  region  of  the  circular  jet,  excited  by  a 
pure  tone.  Based  on  the  frequency  which  gave  maximum  exci- 
tation and  which  was  identical  to  the  measured  vortex- 
shedding  frequency,  the  Strouhal  number  was  also  found  to 
follow  the  relationship  St  = 0.012 /Re  . 

On  dealing  with  the  scalar  mixing  mechanism  in  shear 
flows,  Fiedler  [1974],  described  experiments  conducted  in 
two-dimensional  shear  layers  between  still  ambient  air  and 
a heated  stream  of  air.  Measurements  of  the  cliaracteristics 
of  mean  and  fluctuating  temperatures,  suggested  a predomin- 
ant bulk  conv'jction  mechanism  for  the  temperature  whicl;  was 
explained  by  the  presence  of  large-scale  coherent  structures. 
The  measurements  were  performed  using  conventional  averaging 
and  by  averaging  separately  over  the  turbulent  and  non- 
turbulent  parts  of  the  flow. 

It  was  pointed  out  before  that  considerable  motiva- 
tion in  studying  tlie  mixing  region  of  tvjrbulent  jots  has 
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been  in  its  connection  with  the  noise  generntion  problem. 

The  literature  in  this  topic  is  rather  extensive  and  only 
some  of  the  outstanding  reports  will  be  referred  to  hero. 
Bradshaw  et  al.  [1963]  conducted  a detailed  investigation  in 
a circular  jet  with  an  exit  Mach  number  of  0.3.  They  found 
that  the  flow  in  the  noise-producing  region  is  dominated 
by  a group  of  large  eddies;  these  eddies  in  turn,  appeared 
to  be  the  source  of  the  pressure  fluctuations  near  the  noz- 
zle. They  even  commented  on  the  possibility  of  artificial 
augmentation  of  the  large  eddies  in  order  to  increase  the 
mixing  rate  and  thus  allow  the  reduction  of  noise. 

Mollo-Christensen  and  coworkers,  as  referred  to  by 
Mollo-Christensen  [1967],  have  performed  experiments  in  a 
higher  speed  but  still  subsonic  circular  jet.  Pressure 
crosscorrelation  measurements  suggested  an  ordered  struc- 
ture underneath  the  mixing  layer. 

In  a series  of  publications  (see  for  instance  Kc  and 
Davies  [1971],  Lau  et  al.  [1972],  Fuchs  [1972]),  researchers 
at  the  University  of  Southampton  have  reported  on  studios  in 
regions  adjacent  to  the  mixing  layer  of  a circular  jet.  By 
measuring  correlations  between  pressure  and  velocity  fluctua- 
tions in  both  the  potential  core  and  the  entrainment 
region,  Lau  et  al.  [1972]  inferred  events  occurring  in  the 
adjacent  mixing  layer.  In  fact,  they  postulated  a model  for 
the  structure  of  this  region,  consisting  of  a train  of  dis- 
crete vortices,  evenly  spaced  and  traveling  with  a velocity 
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of  approximately  60%  the  nozzle  exit  velocity  and  a char- 
acteristic dimensionless  frequency  fd/U^  = 0.66. 

As  an  extension  of  the  above  mentioned  investigation, 
Fuchs  [1972]  completed  a survey  of  full-band  and  narrow- 
band  pressure  space-correlations  obtaining  a strong 
coherence  both  laterally  in  a cross-section  and  longitudi- 
nally up  to  eight  diameters  from  the  nozzle  exit. 

Wooldridge  et  al.  [1972]  reported  spectral  and 
correlation  measurements  of  the  axial  velocity  component  in 
the  initial  and  transitional  regions  of  a round  jet  at  vari- 
ous subsonic  regimes.  Their  results  agreed  with  those  ob- 
tained by  Ko  and  Davies  [1971]. 

More  recently,  Lau  and  Fisher  [1975]  conducted  a 
probability  analysis  of  the  time  intervals  between  succesive 
downward  spikes  in  the  longitudinal  velocity  fluctuations 
in  the  region  close  to  the  potential  core  of  the  mixing  lay- 
er. The  highest  probability  corresponded  to  the  vortex 
passing  frequency  found  in  their  previous  investigation. 

Using  a time-domain  averaging  (or  eduction)  technique,  with 
the  spikes  as  triggering  events,  they  determined  educed 
velocity  signals  which  supported  the  proposed  vortex  train 
model . 

Experiments  in  the  wake  behind  a circular  cylinder 
have  been  conducted  by  Bevilaqua  [1973]  and  by  Papailiou 
and  Lykoudis  [1974].  The  vortices  shod  by  a moving  cylinder 
in  still  mercury,  were  found  to  persist  at  high  Reynolds 
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numbers.  An  analysis  of  this  observation  v/as  done  in  con- 
nection to  the  entrainment  problem. 


1 

I 


A comprehensive  review  of  the  flow  around  a circular 
cylinder  has  been  reported  by  Morkovin  [1964].  This  includ- 
ed experimental  observations  and  theoretical  studies  of  var- 
ious aspects  of  the  problem  over  an  extended  Reynolds  number 
range . 


In  summary,  there  is  considerable  evidence  that  the 
mixing  regions  of  turbulent  free  shear  flows  may  be  char- 
acterized by  the  presence  of  relatively  coherent  vortex-like 
structures.  The  passing  frequency  of  these  organized  pat- 
terns when  expressed  through  the  Strouhal  number  fd/U,  takes 
values  of  about  0.6  for  round  jets  and  of  tho  order  of  1 for 
plane  jets  (there  are  fewer  experimental  observations  for 
the  latter,  though) . 

There  are  no  reports  in  the  literature  on  how  the  two 
mixing  layers  of  a plane  jet  interact  when  merging  together, 
nor  on  how  they  may  affect  the  far  field  self-preserving 
region. 

3 . 3 Externally  Disturbed  Plane  Jets 

Besides  being  capable  of  generating  noise,  turbulent 
jets  are  known  to  be  sensitive  to  acoustic  inputs.  Indeed, 
the  latter  effect  when  considering  a flow  v/ith  low  Reynolds 
number  in  tho  verge  to  become  turbulent  has  been  knovm  for  a 
long  time,  as  observed  by  Leconte  [l858] . At  higher 
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Reynolds  nunibers,  the  effects  go  fron  altering  the  process 
of  vortex  forraation,  growth  and  coalescence  in  the  initial 
region,  as  described  by  Rockwell  [1971]  and  Becker  and 
Massaro  [l968]  , to  a direct  interaction  between  the  acoustic 
field  and  the  turbulent  flow. 

A comprehensive  report,  including  their  own  experi- 
ments with  vv^ater  and  air  jets,  was  published  by  Chanaud  and 
Powell  [l962]  . In  a more  recent  investigation,  Simoox  and 
Hoglund  [l97l]  theoretically  analized  the  acoustic  field- 
turbulent  flow  interaction  for  the  plane  jet.  Glass  [1968] 
reported  experiments  performed  in  a supersonic  air  jet  where 
the  sound  waves  generated  in  downstream  positions  affected, 
through  an  acoustic  feedback  process,  the  initial  portions 
of  the  jet,  changing  the  effective  exchange  coefficients  and 
therefore  the  rate  of  jet  spread  and  decay. 

Goldschmidt  and  Kaiser  [l97l]  conducted  experiments 
applying  pure  acoustic  waves  of  various  frequencies  to  the 
self-preserving  region  of  a plane  turbulent  jet.  Their  con- 
cluding results  were  that  self-preservation  of  the  flow  is 
maintained,  and  that  at  specific  frequencies,  the  spreading 
rate  is  increased  (with  an  accompanying  faster  decay  of  the 
centerline  velocity)  and  the  geometric  virtual  origin  is 
shifted. 

Considerable  interest  in  the  above  mentioned  aspects 
of  turbulent  jets  has  arisen  from  the  current  technology  on 
fluidic  control  devices.  Particularly,  in  respect  to  the 


dynamic  response  of  fluidic  amplifiers,  the  loss  of  perform- 
ance at  high  frequencies  is  believed  to  be  related  to 
changes  of  the  turbulent  jot  characteristics  when  subjected 
to  periodic  initial  disturbances.  This  was  pointed  out  by 
Stiffler  [1971]  in  his  study  of  the  dynamic  characteristics 
of  the  mean  velocity  profile  of  a plane  jet. 

Binder  and  Favre-Marinet  [l972]  reported  experiments 
done  on  a plane  jot  subjected  to  relatively  large  amplitude 
transverse  blowing  pulses  in  the  nozzle  near  region.  The 
pulses  were  generated  with  a cylindrical  rotating  pressure 
chamber,  provided  with  a blowing  slot  on  each  side  of  the 
nozzle,  and  were  alternatively  periodic.  A higher  rate  of 
spreading  and  higher  mean  velocity  decay  at  the  centerline 
were  found,  when  compared  with  the  undisturbed  jet. 

It  is  seen  therefore,  that  turbulent  jets  are  sensi- 
tive to  a considerable  extent,  to  external  disturbances, 
either  in  the  form  of  pure  acoustic  tone  or  as  localized 
gross  blowing  inputs.  The  initial  region,  of  the  jet  is  par- 
ticularly affected  resulting  perhaps,  in  changes  of  any 
possible  relationship  between  the  coherent  structures  in 
this  region  and  the  flapping  behavior  far  dov/nstream  of  the 
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CHAPTER  4 

EXPERIMENTAL  SETUP  AND  PROCEDURES 

4 . 1 Experimental  Setup 

All  the  experiments  and  measurements  pertaining  to 
this  work  were  done  in  one  of  the  existing  turbulent  plane 
jet  setups  at  the  Ray  W.  Herrick  Laboratories  of  Purdue 
University.  In  the  next  paragraphs  a description  of  the  ex- 
perimental setup,  which  was  essentially  the  same  used  by 
Ott  [1972]  and  Young  [1973],  is  followed  by  a detailed  ac- 
count of  the  instrumentation  and  measurement  procedures  used 
in  this  investigation. 

The  facility  consisted  of  a vertical  rectangular 
nozzle  of  0.635  x 30.48  cm  (aspect  ratio  4 8 ) supplied  with  air 
by  a 1/4  HP  squirrel-cage  blower.  A flexible  duct  connector 
between  the  blower  outlet  and  the  plenum  chamber  was  employed 
in  order  to  avoid  the  transmission  of  vibrations  from  the 
blower.  The  plenum  chamber  of  dimensions  10.25  cm  width, 

30.48  cm  height  and  67.5  cm  long,  was  provided  with  a series 
of  honeycombs  and  screens  to  ensure  uniform  flow.  The  test 
section  downstream  of  the  nozzle  exit  was  confined  between 
two  90.2  X 122  cm  plywood  horizontal  plates  (Figure  4.1). 

In  this  way,  a turbulent  jot  with  Reynolds  number  based  on 
the  slot  width  (0.635  cm)  and  exit  velocity  (2 . 4 3 x 10^  cm/sec) 
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of  approximately  10  , was  obtained. 

The  existing  two-dimensional  traversing  mechanism, 
with  range  152.4  cm  in  the  longitudinal  x-direction  and 
76.2  cm  in  the  lateral  y-direction,  was  modified  by  instal- 


ling two  10  turns  precision  potentiometers  for  lateral  in- 
dependent positioning  of  the  probes.  (Figure  4.2).  Once 
the  potentiometers  were  calibrated,  positioning  of  probes 


within  ± 0.01  cm  was  obtained. 


In  order  to  measure  mean  velocity  profiles  and  to 
check  self-preservation  of  the  flow,  a pair  of  Pitot-tubes 
(0.1  cm  internal  diameter)  were  manufactured.  The  Pitot- 


tubes  were  used  together  with  a pressure  transducer  and  a 
X-Y  plotter,  the  former  previously  calibrated  against  a 
micromanometer  (Figures  4.3  and  4.4). 


Traverses  were  done  at  various  longitudinal  posi- 
tions. The  resulting  mean  velocity  profiles  appear  in 
Appendix  A,  where  the  self-preserving  characteristics  of  the 
flow  are  also  presented. 

Two  DISA  55A01  constant-temperature  anemometers  with 
DISA  55A22  and  DISA  55F11  probe  supports  and  "home-made" 
plated  tungsten  sensing  elements  (0.0002  cm  diameter),  were 
used  to  obtain  the  velocity  signals  at  any  two  points  in  the 
flow.  In  computing  normalized  crosscorrelations  between 
these  signals,  calibration  of  the  hot-wires  v/as  excluded. 


It  sufficed  to  measure  the  root-mean-square  value  of  each 
signal  and  divide  as  noted  in  equation  (2.17).  A true  root- 
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moan-squaro  voltmeter  R&K  215  was  used  for  this  purpose. 

Preliminary  crosscorrelation  function  measurements 
were  done  using  a SAICOR  Correlator  and  Probability  Analyzer, 
model  SAI-42.  This  digital  on-line  processing  instrument, 
provides  a correlation  computation  at  100  incremental  lag 
points  simultaneously.  By  using  the  resume  capability  of 
the  instrument,  an  average  of  various  correlation  samples 
can  be  determined,  although  the  process  is  slow  and  tedious. 

In  order  to  overcome  the  above  disadvantage  and  to 
have  additional  computational  features,  the  use  of  a Hewlett- 
Packard  5452A-2114B  Fourier  Analyzer  System,  was  preferred. 

It  consists  of  the  following  elements:  an  analog-to-digital 
converter  with  a wide  range  of  sampling  characteristics;  an 
8K  words  memory  minicomputer  which  includes  the  Fast  Fourier 
Transform  algorithm  among  its  software;  and  ^ display,  a tele- 
printer and  a fast  paper  tape  puncher  units.  Depending  on 
the  data  block  size  used,  a more  complete  picture  of  the 
correlation  function  could  be  obtained,  that  is,  including 
more  lag  points.  Also  by  programjning  the  system,  any  number 
of  samples  could  be  taken  for  a particular  situation  and 
therefore  a better  estimate  of  the  correlation  function  was 
obtained.  Detailed  block  diagrams  and  listings  of  the  pro- 
grams used  in  this  investigation  are  presented  in  Appendix 
B. 

Finally,  auxiliary  equipment  like  two  SKL302  low- 
pass  filters  to  avoid  aliasing  in  the  analog-to-digital 


conversion,  a KISOO  X-Y  Plotter,  a general  ourpose  oscillo- 
scope and  digital  voltmeters,  \;ere  also  used.  (Figures  4,5, 
4.6,  4.7). 

4 . 2 Experimental  Procedures 

The  measurement  procedure  v/as  essentially  the  same 
in  all  experiments.  After  having  carefully  positioi^ed  and 
aligned  the  probes,  a warm  up  period  (for  the  instrumenta- 
tion) of  at  least  one  hour  was  allowed.  At  tl\e  beginning  of 
each  run,  a quick  average  correlation  based  on  10  or  20  sam- 
ples was  determined  with  the  Fourier  analyzer.  In  this  way, 
proper  scaling  and  sampling  parameters  for  the  annlog-to- 
digital  converter  were  determined  (the  cut-off  frequencies 
or  the  low-pass  filters  were  set  to  just  avoid  aliasing. 
Generally  1 kHz,  but  no  less  than  GOO  Hz,  was  used  for  these 
frequencies).  After  this  preliminary  operation  the  analyzer 
was  reprogrammed  and  up  to  four  estimates  of  the  correlation 
function  were  determined  for  each  run  at  each  point  of 
interest.  Each  estimate  included  250  samples  and  took  about 
20  minutes  to  be  obtained. 

Four  kinds  of  velocity  correlation  measurements  were 
taken.  These  are  labelled  Type  I,  II,  III  and  IV,  respec- 
tively. Figures  4.8,  4.9,  4.10  and  4.11  show  schematically 
how  the  probes  were  positioned  for  each  measurement  typo, 
respectively.  Table  4.1  gives  the  various  values  of  the  x 
and  y coordinates  used  in  positioning  the  probes  for  each 
type  of  measurement. 
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MoasuromontG  Typo  I (Fiyuro  4.8)  consisted  of  cross- 
correlations  between  x-componont  velocity  signals  obtained 
from  probes  set  at  two  points  symmetrically  apart  with 
respect  to  the  centerline  of  the  jet.  These  measurements 
were  performed  at  x/d  stations  of  10,  20,  30,  40,  60,  80  and 
100;  the  probes  were  positioned  at  y/b  values  of  0.5,  0.75, 
1.0  and  1.25  approximately,  for  each  x/d  station. 

Corresponding  to  the  Type  II  runs  (Figure  4.9), 
crosscorrelations  were  measured  by  setting  tlie  probes  at 
points  with  different  opposite  lateral  coordinate  y/b.  The 
Type  II  runs  were  limited  to  x/d  = 40.  However,  a wide 
number  of  combinations  for  the  lateral  coordinates  y^^/b  and 
y^/b  wore  included. 

The  third  typo  of  measurements.  Type  III  (Figure 
4.10) , was  obtained  by  holding  one  probe  at  a given  x/d  sta- 
tion and  placiiig  the  second  probe  at  the  opposite  side  of  the 
jet  but  at  a different  x/d  position.  The  longitudinal  sep- 
aration between  probes  Ax,  was  varied  from  0 to  5.1  cm  in 
steps  of  1.27  cm.  Three  x/d  stations  (20,  40  and  GO)  for 
the  fixed  probe  wore  investigated. 

Finally,  the  measurement  procedure  of  tl\c  Type  IV 
(Figure  4.11)  was  iJifferent  from  those  .ilready  doscri)>ed . For 
this  case  the  autocorrelation  function  of  the  y-componont 
velocity  at  the  centerline  of  the  jot,  was  measured  for  x/d 
positions  of  20,  30,  40,  60  and  80.  Instead  of  two  normal 
hot-wires,  an  x-wire  typo  DISA  55A38  was  omiJloyod  for  those 
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measurements.  Detailed  analysis  on  how  the  autocorrelation 
was  estimated  for  this  case  is  yiven  in  Appendix  B,  toijether 
with  the  program  used  to  process  the  data  with  tlie  Fourier 
analyzer. 

Measurements  wore  also  taken  at  the  mouth  of  the  jot 
to  attempt  to  determine  any  instabilities  in  the  blower 
feeding  the  nozzle.  Preliminary  and  liiritcd  observati  tins 
were  also  made  in  the  presence  of  a low  t'roepiency  acoustic 
disturbance.  These  are  noted  in  Chapter  7. 

As  a last  series  of  experiments,  modification  of  the 
electrical  motor  and  pulleys  driving  the  blower,  allowed 
runs  at  Reynolds  numlscrs  of  7900  and  15100.  These  included 
measurement  of  crosscorrelations  at  x/d  stations  of  40,  60 
and  80,  with  the  probes  synunetrically  positioned  with  re- 
spect to  the  jet  centerline  (Type  I measurements) . 

The  objetive  in  all  the  measurements  (Types  I to  IV) 
was  to  first  determine  whether  indeed  there  is  a flapping 
motion  of  the  jot  and  secondly  to  attomj^t  to  diagnose  its 
origin. 
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CHAPTER  5 

EXPERIMENTAL  RESULTS 

In  this  chapter  the  experimental  results  are  pre- 
sented. In  order  to  have  a better  depiction  of  the  flapping 
behavior  first  some  general  properties  are  reported  and  then 
various  characteristics  (as  resulted  from  the  experiments) 
are  separately  considered.  The  various  types  of  measure- 
ments, tabulated  on  Table  4.1,  are  reported  in  the  sections 
that  follow  as  needed  rather  than  chronologically.  Generally 
the  type  I measurements  were  taken  to  determine  flapping 
frequency  and  amplitude.  The  type  II  measurements  were 
taken  to  assure  that  the  motion  was  transversely  coherent. 

The  absence  of  "weaving"  of  the  flow  past  a fixed  vortex 
street  type  structure  was  noted  through  type  III  measure- 
ments. They  indeed  confirmed  the  uniform  sideways  motion  of 
the  entire  flow  field  and  gave  an  approximate  wave  velocity 
(or  convective  velocity) . The  type  IV  measurements  were 
taken  to  insure  that  the  sideways  motion  was  also  noted  at 
the  jet  centerline. 

Additional  measurements  at  different  Reynolds  num- 
bers are  reported  in  the  last  part  of  this  chapter. 


5.1  Charactoristics  of  tho  IMappincj  r^otion 

All  nKvu^urcd  crossoorielation  functions  between 

.K-components  of  the  velocities  at  two  points  on  opposite 
sides  of  the  centerline  of  the  jet,  are  persistent  and  re- 
peatable. These  correlation  estimates  present  some  general 
features  which  will  bo  examined  in  the  following  paragraphs. 

The  two  most  important  characteristics,  indicating 
the  existence  of  an  oscillating  (in  the  average)  lateral 
motion  or  flapping  of  the  jot,  are:  an  invariably  negative 
value  of  the  correlation  function  at  cero  time  delay  and  an 
alternate  succession  of  positive  and  negative  values  as  the 
time  delay  between  the  two  sigiials  increases.  Examples  of 
those  correlations  (as  obtained  directly  from  the  Fourier 
analyzer  tlirougli  an  x-y  plotter),  are  given  in  Figures  5.1 
to  5.*).  The  correla tion  functions  shown  in  Figures  5.1  and 

5.2  were  obtained  with  tlie  i^ot-wire  probes  positioned  sym- 
metrically apart  with  respect  to  the  jet  centerline  at 
.\/d  = 40,  y/b  = t 1.0  and  x/d  = 100,  y/b  = + 1.0,  respec- 
tively (Type  I).  Figures  5.3  and  5.4  sliow  correlations 
taken  at  x/d  = 40,  witli  the  probes  asymmetrically  positioned 
(Type  II):  Figure  5.3  corresponds  to  y^/b  = 0.75,  Y 2^^  ~ 

- 0.5  (Note  tliat  for  t > 0.128  sec  in  Figures  5.1  ami  5.2 
and  for  r > 0.064  sec  in  Figures  5.3  and  5.4,  that  is,  the 
second  half  of  all  records,  the  correlation  function  is  not 
significant.  This  is  ikie  to  the  wrap-around  error  inherent 
in  the  analyzer  piU'formance , as  explained  in  Appemlix  13), 


Crosscorrelation  Function  R(T)  [volts 


Crosscorrelation  Function  R(T)  [volts 


Moasuroments  Typo  II 


Time  Delay  T 


An  Example  of  Crosscorrolation  Function 

y /b  = 1.25,  y,/b  = - 


Figure  5.-1 
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The  correlations  exhibit  a certain  periodicity. 
This  is  related  to  and  gives  a measure  of  the  flapping  fre- 
quency. (Simultaneously,  Uie  magnitude  of  the  correlation 
at  I = 0 would  give  a measure  of  the  amplitude  of  the 
oscillatory-like  flapping  motion) . The  times  between  local 
maxima  or  local  minima  (see  for  instance  Figure  5.1)  are 
generally  not  equal.  However,  the  time  to  the  first  maxi- 
mum is  clearly  measurable  and  it  suggests  the  definition  of 
a frequency  of  flapping,  f^^,  as  illustrated  schematically  in 
Figure  5.5.  On  the  other  hand,  an  average  time  between 
local  maxima  and  minima,  beyond  the  first  local  maximum,  can 
be  found.  With  it  a second  frequency  f^,  as  also  illustra- 
ted in  Figure  5.5,  can  be  determined.  Those  two  frequencies 
are  essentially  the  sam.e , although  within  some  experimental 
scatter,  as  noted  in  the  beginning  of  the  next  section. 

5 . 2 Characteristics  of  the  Flapping  Frequency  along  the 
Longitudinal  and  Lateral  Coordinates 

Type  I Measurements 

A complete  survey  of  crosscorrelation  functions  was 
done  after  setting  both  probes  on  opposite  sides  of  the  jot 
centerline  (measurements  Type  I).  The  resulting  frequencies 
of  flapping  arc  plotted  as  a Strouhal  number,  fd/U^.  Fig- 
ures 5.6  and  5.7  present  the  results  for  f^  and  respec- 

tively. The  flapping  frequency  (given  by  f^^  and  f2)  is 
noted  to  decrease  as  x/d  increases. 


Strouhal  Number 


Dimonsioiilous  I,ciinjitud  inal  Coord iiiaLo  x/d 


Kitjuro  5.7  Strouhal  Numbor  f^d/U^ 


x/d 


VS  . 
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On  the  other  hand,  the  flapping  frequency  does  not 
seem  to  depend  on  the  lateral  coordinate.  This  is  noted 
both  by  the  collapse  of  the  data  points  on  Figures  5.6  and 
5.7  and  by  the  plot  in  Figure  5.8.  This  is  so  at  least  for 
x/d  > 30  where  the  jet  becomes  self-preserving  in  terms  of 
mean  velocity  and  turbulence  intensity  (Bradbury  [1965], 
Gutmark.  and  Wygnanski  [1976]).  Tlie  data  points  correspond- 
ing to  x/d  values  of  10  and  20  show  a large  amount  of  scat- 
ter. The  reason  for  this  could  be  error  in  positioning  the 
probes  at  these  small  x/d  stations  and  insufficient  time 
resolution  in  the  correlations  computed  by  the  analyzer,  as 
compared  with  larger  x/d  points  at  lower  velocities. 

Several  scaling  parameters  for  f (both  f^  and  ^2^ 
wore  considered.  The  intent  was  to  exhibit  any  possible 
similarity  and  self-preservation  for  the  flapping  frequency. 
The  various  jet  coordinates,  length  and  velocity  scales  (x, 
y,  d,  b,  U,  U ) were  combined  in  this  search.  The  best  re- 
suits  are  shown  in  Figures  5.9  to  5.15. 

Figure  5.9  shows  the  variation  of  — q—  against  y/d 
for  each  .x/d,  where  U is  the  local  mean  velocity  at  each 
pair  of  points  of  coordinates  x/d,  i y/d  (where  the  cross- 
correlation was  measured).  Figure  5.10  presents  a collapse 
of  all  curves  in  Figure  5.9,  using  y/b  instead  of  y/d  for 
the  abscissa. 

The  same  data,  but  in  somewhat  different  form,  is 
presented  in  Figure  5.11  when  the  parameter  is  plotted 
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Figure  5.12  Dimensionless  Flapping  Frequency  f,y/U 
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acjain-st  (y/b)  , using  a logarithmic  scale  for  the  ordinate. 
Tlio  data  points  follow  approximately  a straight  line,  that 
is : 


In 


fX/u 

y/b 


(5.1) 


or 


^ = 0.11  ^ exp  0.72  (^) 


(5.2) 


Figure  5.12  relates  fy/U  (where  U is  the  center- 

m m 

line  mean  velocity  at  a given  x/d)  to  y/b.  A fit  of  the 
type: 


0.11  ^ 

b 


(5.3) 


appears  satisfactory. 

A.ctually  both  equations  (5.2)  and  (5.3)  represent 
the  same  trend.  In  effect,  if  Reichardt's  solution  for  the 
self-preserving  mean  velocity  profile  is  recalled  (equation 
(2.4) , Chapter  2) : 


^ = exp  (-  (0.8  326  ^-)^J 


^ = exp  [-  k (^d  ^1 

m 

it  may  be  combined  with  (5.2)  to  give  (5.3). 


(2.4) 


(5.4) 


Moreover,  equation  (5.3)  may  also  bo  written  as: 


j.j  nil  oi  crosscorreiacion  Function  roi*  x/d  = 40, 

Vj^/b  = 0.75,  y^/b  = - 0.25 


6 2 


= 0.11  (5.5) 

m 

which  is  the  result  sliown  in  Fiyures  5.13,  5.14  and  5.15. 

The  flapping  frequency  then,  attains  approximate 
self-preservation  if  scaled  with  the  jot  half-width  b and 
the  mean  velocity  at  tlie  centerline. 

Using  this  result,  a replotting  of  the  correlation 
functions  was  done  througli  a computer  Calcomp  system.  The 
abscissa  i was  modified  to  a dimensionless  time  delay 
I U^/b  iind  the  ordinate  was  divided  by  the  product  of  the 
r.m.s.  values  of  the  signals  (to  have  a normalized  crosscor- 
relation function).  Figures  5.16  to  5.19  show  tlicso  results 
for  y/b  values  of  0.5,  0.75,  1,  and  1.25  respectively,  each 
one  including  x/d  stations  larger  than  30  (tl\c  plots  for 
x/d  = 10  and  20  were  purposely  omitted.  They  did  not  sliow 
the  siime  trend)  . 

It  is  interesting  to  note  that  not  only  does  hlic 
flapping  frequency  exhibit  approximate  self-preservation  and 
similarity  for  x/d  > 30,  but  so  does  the  magnitude  of  the 
crosscorrelation  function  (at  least  for  the  first  half- 
period) . As  an  example,  the  values  of  R(0),  where: 


u(y  , t)  u(-y  , t) 

R(0)  = 

/ — / — 

/ u (y,t) / u“ (-y,t) 


(5.6) 


are  plotted  in  Figure  5.20. 
ship  of  tlie  form: 


A trend  is  noteil  and  a relation 


x/d  y/b 


Figure  5.17  Normalized  Crosscorrelation  Functions,  y/b 
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U(0)  (1-k^  y/1^)  ‘■■•xp  (-k^  y/b)  (5.7) 

seems  to  be  satisfied.  Self-preservation,  however,  is  far 
from  perfect.  In  this  respect  it  should  be  remembered  tliat 
the  more  complicated  a statistici\l  quantity  is,  the  longer 
it  takes  to  reac)\  self-preservation  (Tennekes  and  Lumlcy 
[1972]).  The  best  fits  for  the  coefficients  in  equation  5.7 
appear  to  be  k^  = 2.5,  k2  - 2.9.  'I'liose  are  not  shown  on 
Figure  5.20.  It  does  make  sense  to  expect  that  for  sm.ill 
y/b  the  typical  lateral  correlation  function  is  obtained 
(Hinze  [1959]),  whereas  for  y/b  ''  0.35  tlie  neijativ’C  values 
attributed  to  flapping  become  predominant. 

A detailed  survey  of  correlation  fr.nctions  v;as  per- 
formed with  both  probes  at  x/d  = 40  but  at  equal  |y/b|  loca- 
tions. This  was  done  in  order  to  have  a better  picture  of 
the  flapping  beliavior  across  a given  longitudinal  coordinate 

and  extend  the  results  of  Figure  5.20.  Figures  5.21  a,  b,  c, 
c,  e and  f show  a sequence  of  such  functions  for  and  increasing 

separation  distance  between  the  two  symme tr i ea  1 ly  se.t  probes. 

Tlio  curves  corresponding  to  ± y/b  < 0.35  shown  in  Figure  5.21 

a,  b,  and  c,  presented  a positive  value  for  h(0) . Note  however, 

that  R(0)  decreases  as  the  separation  distance  between  [jrobes 

increases  (this  is  an  expected  result  as  commented  later). 

Moreover,  at  t y/b  = 0.35  the  value  of  K(0)  is  smaller  that 

the  values  of  R(t)  corresponding  to  slightly  positive  and 

negative  time  delays  (Figure  5.21c).  At  ' y/b  = 0.4,  R(0)is 

aproximately  zero  aiul  the  correlation  function  taf.es  a kind 


rosscorrelation  Functi 
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of  transitional  shaiJO  (Figure  5.7.  Id).  Finally,  for 
1 y/b  > 0.‘15,  the  correlations  characterizing  the  flapping 
beliavior  were  obtained  (Figures  5.21  e and  f)  . The  values 
of  R(0)  in  Figures  5.21  correspond  to  the  values  in  the  plot 
of  Figure  5.20. 

lYpc  11  Measurements 

It  was  mentioned  before  that  the  fla[jpinf/  fr<?quoncy 
did  not  e.vhibit  a preferred  dependeiice  on  tlie  lateral  co- 
ordinate (Figure  5.8).  This  was  one  of  the  results  of  the 
Type  I measurements  (those  performed  with  tlie  probes  syrruiiet- 
rically  positioned  with  respect  to  tlic  jet  centerline)  . 

A question  may  still  arise  as  to  whether  the  motion, 
interpreted  as  flapping,  is  instead  due  to  a "weaving"  of 
the  mean  flow  through  an  impressed  vorto:  type  structure. 

To  determine  that  possibility  Type  II  measurements  were 
talcen  with  the  sensors  at  different  lateral  locations  (sam- 
ple plots  are  shown  in  Figures  5.3  and  5.4) . The  presence 
of  any  such  structures  v.'ould  have  to  give  further  changes  in 
the  sign  of  the  correlation  at  x 0.  This  v/as  not  noted. 
Fur tl'.ermore  no  dependence  of  tlie  flapjjing  frequency  on  lat- 
eral location  was  observed  eitlier. 

Tlie  laisults  of  measurcHl  flapiiing  freriucncy  are  com- 
bined in  Figure  5.22  (compare  willi  Figures  5.13,  5.14  and 

5.15).  The  average  value  of  the  oarameter  fb/U  is  constant 

‘ m 

and  compares  well  with  that  obtained  from  the  Type  I 


measurements . 
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It  seems  therefore,  that  the  frequency  or  flapping 
when  expressed  as  fb/U^^,  is  a single  value  characteristic 
for  each  x/d  station,  independently  of  the  points  considered 
in  measuring  the  crosscorrelation,  for  |y/bl>_  0.5. 

The  amplitude  of  the  correlation  function  at  zero  time 
delay  R(0) , is  however,  a function  of  y/b  as  shown  in  the 
sequence  of  correlations  of  Figure  5.21.  This  is  an  expect- 
ed result,  since  R(0)  as  defined  and  obtained  by  varying  the 
separation  distance  between  probes,  is  equivalent  to  the 
space-correlation ; 

R(y,Ay)  = uTyl-'.Ay)  u(y-TML  (5.8) 

y^u^(y+SAy)  y^u^(y-'5Ay) 

for  y = 0.  In  equation  (5.8),  y = S (y2^+y2)  ^nd  Ay  is  the 
|i  separation  between  probes. 

Further  analysis  was  done  of  the  Type  I and  Type  II 
il  data  at  x/d  = 40.  The  objective  was  to  obtain  space- 

correlations  in  a form  similar  to  equation  (5.8).  The  re- 
sults are  presented  in  Figure  5.23.  The  lateral  coordinate 
of  the  middle-point  between  probes  (y  = S (y^+y2^  ) 
separation  distance  Ay  = y2^“y2'  nondimensionalized  by 

the  jet  half-width  b.  Mote  that  the  space-correlation  func- 
tion never  attains  a negative  value  for  y/b  1 • 0 / 
approaches  zero  asymptotically  as  the  probe  separation  Ay/b 


I 


\ 

i 

i 


increases . 
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The  same  curves  are  presented  in  Figure  5.24  with  a 
smaller  scale  in  order  to  compare  them  with  similar  re- 
sults reported  by  Gutmark  and  Wygnanski  [1976].  Further 
comments  on  this  comparison  are  presented  in  Chapter  6. 

5 . 3 Amplitude  of  Flapping 

An  estimate  of  the  amplitude  of  flapping  can  be  done 
by  assuming  a sinusoidal  bodi]y  displacement  of  the  moan 
velocity  profile  of  the  jet  from  its  mean  location. 

In  effect,  the  instantaneous  displacement  of  the 
mean  velocity  profile  may  bo  assumed  as: 

e = sin  2TTft  (5.9) 

where  f is  the  frequency  of  flapping  and  is  the  ampli- 
tude to  be  estimated. 

If  e is  the  instantaneous  a.c.  voltage  of  the  hot- 
wire anemometers,  tlion  for  an  essential  linear  response,  the 
correlation  function  can  be  expressed  as: 


_ e(x,y,t)  e(x,-y,tt-f) 
/^e^(x,y,t)  y^e^  (x,-y  , tFT ) 


(5.10) 


Furthermore,  for  sufficiently  small  values  of  r,  the  nonran- 
dom component  of  e can  be  approximated  by: 


(5.11) 


where  E is  the  d.c.  voltage  output  of  the  anemometers. 
Substituting  (5.9)  in  (5.11)  and  (5.11)  in  (5.10),  one 


Gutmark  and  Wygnanski  [1976 


Figure  5.24  .Normalized  Lateral  Correlation  R{y,Ly)  vs.  Ay/b 


obtains : 


R (t  ) = 


(5.12) 


On  the  other  hand,  an  empirical  cooling  law  for  the 
hot-wires,  may  be  expressed  as: 


R (R  -R.)  ^2^ 

w w f 


2 Si 

E = A + BU^ 


(5.13) 


(5.14] 


where  A = k^R^  (R^^-R^ ) , B = k2R^^(R^^  -Rf)  , and  k2  are  em- 
oirical  constants  and  R and  R^  are  the  operating  and  cold 
resistance  of  the  hot-wire,  respectively. 

Differentiating  (5.14)  with  respect  to  U: 


2e|£  = SBU-’ 


3U  4'  E 

combining  (5.14)  and  (5.15),  one  gets: 


(5.15) 


M = -1- 

3U  4U  E 


(5.16) 


Thus,  writing  substituting  in  (5.12)  after 

using  (5.16)  : 
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R(t)  = - 


r*  2 ~i  ^ 

1 E -A  nu 

1 TU  “iT-  Sn  2 


:5.i7) 


At  T = 0,  the  last  expression  reduces  to: 

r 1 E“-A  Dul  2 


R(0)  = I 


■> 

I E“-A  3U 

^ ^ 3y_ 


from  which,  a relationship  is  obtained  for  c 


= [2R(0) 

^ 1 E^-A  au 


4U  E ay 


(5.18) 


Furthermore,  expressing  c , U and  - — in  dimensionless 

^ m ay 

form  by  using  the  half-width  b and  the  mean  velocity  at  the 
centerline  of  the  jet,  (5.18)  reduces  to: 


[2R(0) P 

U E^-A  a(u/u  ) 
m m 

4U  E a(y/b) 


(5.19) 


Finally,  assuming  a mean  velocity  profile  as  given 
by  Reichardt  (equation  2.4,  Chapter  2): 


= exp  -(0.8326 
U b 

m 


(2.4) 


and  taking  derivatives,  one  obtains: 


U ^ ^ 

= _ 1.39 


(5.20) 
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Expression  (5.20)  substituted  in  (5.19)  (jives; 


4.07 


'[R(0)] 
E b 


(5.21) 


Limited  data  was  ta):en  which  included  d.c.  voltage 
(E)  values  for  the  anemometers.  An  approximate  value  for 
the  constant  A was  estimated,  by  extrapolating  to  u’/^  = 0, 
the  straight  line  represented  by  (5.14),  This  was  obtained 
from  measurements  of  the  type  I at  x/d  of  80,  60  and  40. 

The  approximate  Cjjj/b  values  calculated  through 
equation  (5.21)  are  tabulated  in  Table  5.1.  As  noted,  the 
amplitude  of  flapping  is  in  the  order  of  15%  to  23%  of  the 
jet  half-width,  and  does  not  appear  to  depend  on  x/d  nor  on 
y/b. 


Table  5.1  Amplitude  of  Flapping 


x/d 

y/b 

80 

1.13 

i 

0.18 

0.85 

0.21 

0.56 

0.18 

60 

1.07 

0.17 

0.80 

0.19 

0.53 

0.23 

40 

1 

1.03 

0.15 

0.77 

0.16 

80 


5 . 4 Convection  Characteristics  of  the  Flapping  Behavior 

A series  of  correlation  measurements  were  done  with 
the  hot-wire  probes  at  different  longitudinal  coordinates  on 
opposite  sides  of  the  jet  (Type  III  measurements).  The 
original  intent  was  to  further  investigate  the  uniform 
flapping  motion  of  the  jet.  The  data  also  gives  insight  on 
how  the  flapping  behavior  is  convected  downstream  of  a given 
x/d  position.  One  probe  was  set  at  a fixed  x/d  station  and 
the  second  was  successively  positioned  at  steps  of  2d  in  the 
downstream  direction.  The  corresponding  crosscorrelation 
functions  were  determined  according  to  equation  (5.22)  be- 
low: 


r(t)  = u(x/d,y/b,  t)  u ((x/d ) +2n , -y/b , t+xT 


(5.22) 


/ (x/d,y/b,t)  / X?  ( (x/d)+2n,-y/b,t+T) 
withn=  0,  1,  2,  .../  5 


Four  x/d  stations  (10,  20,  30  and  40)  were  surveyed, 
keeping  both  probes  at  the  same  lateral  coordinate  ± y/b  = 1. 
Three  more  (20,  30  and  40)  were  done  with  ± y/b  = 0.65. 

Each  x/d  station  included  downstream  separations  for  the 
second  probe,  up  to  10d(n  = 5).  Figures  5.25  a,  b,  c,  d,  e, 
f and  5.26  a,  b,  c,  d,  e,  f show  examples  of  sequences  of 
such  correlation  functions  for  increasing  downstream  separ- 
ation. Figure  5.25  corresponds  to  x/d  = 10,  ± y/b  = 1.0, 
and  Figure  5.26  to  x/d  = 20  + y/b  = 0.65.  As  noted,  the 
time  delay  corresponding  to  the  negative  probe  (pointed  with 
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‘^rrow)  ivs  tl\c  sopm'iitiofv  bot.wcct\  probes  iiicrocisGs. 

effectively  indicates  thiit  the  fl.\i>pvi\cj  bohcivior  i\s 
detected  at  a given  longitudinal  position  can  be  sensed  at 
a later  time  at  a downstream  position. 

An  indication  of  the  velocity  of  convection  of  the 
flapping  behavior  can  be  obtained  by  plotting  the  time  delay 
shift  in  the  corielation  functions  against  t)io  corresponding 
longitudinal  separation  between  probes.  Fitjures  5.2  7 to 
5.33  show  this  result  for  each  case  investigated. 

Based  on  the  separation  between  probes  Ax  and  the 
time  delay  shift  At,  a velocity  of  convection  of  the  flap- 
ping behavior  can  be  defined  as  follows: 


U 


cf 


Ax 

At 


Ax-*-0 


(5.23) 


This  can  be  roughly  obtained  by  extrapolating  the 
curves  in  Figures  5.27  to  5.33  to  zero  separation  distance 
and  computing  tlie  ratio  in  exiiression  (5.23)  above.  The 
corresponding  flapping  convective  velocities  divided  by  the 
corresponding  local  mean  velocities  (d^,j./U)  , and  by  the 
centerline  mean  velocities  arc  talnilated  in  Table 

5.2. 

Also  indicated  are  the  convective  velocities  (of  the 
turbulent  structure)  measured  by  Young  in  the  same 

jet  setup.  These  velocities  were  obtained  by  measuring 
space-time  correlations  between  longitudinal  velocity  com- 
ponents at  two  points  on  the  same  side  of  the  jet.  The 
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Figure  5.29  Time  Delay  Shift  At  vs.  Ax/d,  y/b  = 1, 
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Table  5.2  ConvocLivo  Velocities  of  I'lappiny 
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x/d 

y/b 

"ctA 

U ..  / U 
c t / m 

* 

10 

1 

1.16 

0.58 

20 

1 

1.28 

0.64 

1.68 

30 

1 

1.15 

0.  57 

1.43 

40 

1 

1.23 

0.61 

1 . 58 

20 

0.65 

0.99 

0.74 

1 . 18 

30 

0.65 

0.94 

0.70 

1.  1 

40 

0.65 

1.01 

0.75 

1.18 

★ 

Moasuremotits  by  Young  (197J]. 


direction.  As  noted  in  Table  5.2  tlio  flapping  convective 
velocities  resulted  about  30'b  smaller  than  the  convective 
velocities  measured  by  Young  for  y/b  = 1,  and  about  20%  for 
y/b  = 0.65. 

5.5  Autocorrelation  of  the  Tiateral  c’oniponcnt  Velocity  at 
the  Centerli  n e of  tl^e  Jet 

All  the  results  presented  in  the  previous  sections 
of  this  chapter,  were  obtaii\cd  for  longitudinal  velocity 
components,  easily  measured  with  normal  wires.  In  the  pre- 
sence of  flapping  the  lateral  velocity  components  should 
show  similar  behavior.  To  obtain  such  u;casurements , two 
x-wire  probes  sliould  be  usctl.  Althou<jh  the  e<]uipmont  was 
available,  the  use  of  four  hot-wire  aiK'mometers  and 
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corresponding  signal  processing  equipment  would  have  led  to 
considerably  larger  averaging  times.  As  a result  only  a 
spot-check  was  made  with  a single  x-wire  probe. 

Autocorrelation  functions  of  the  y-component  of  the 
velocity  were  determined  along  the  jet  centerline,  at  x/d 
positions  of  20,  30,  40,  60  and  80.  (Type  IV  measurements). 

The  data  should  obviously  give  the  normal  autocorrelation 
function  for  small  delay  times,  with  a maximum  (equal  to 
the  mean  square  value  of  the  signal)  at  x = 0.  However,  for 
large  delay  times,  the  flapping  motion  (if  at  all  present) 
should  be  accompanied  by  alternatively  positive  and  nega-  ' 
tive  values  of  the  autocorrelation.  The  results,  with  a 
single  x-wire  suffer  in  signal  resolution  but  do  indeed  show 
the  expected  trends.  Figures  5.34,  5.35  and  5.36,  show  ex- 
amples of  such  functions  as  obtained  directly  from  the 
Fourier  analyzer  through  an  x-y  plotter,  for  x/d  = 20,  40 
and  80  respectively. 

As  noted,  the  autocorrelation  functions  have  their 
maximum  value  at  zero  time  delay.  On  the  other  hand,  the 
curves  are  characterized  by  a pseudo-oscillatory  nature  as 
the  time  delay  increases,  consistent  with  the  previously  de- 
scribed crosscorrelation  functions  of  the  x-component  vel- 
ocity. The  dimensionless  flapping  frequency,  fd/U^,  can  be 
determined  as  before  and  characterized  by  both  fj^  and 
The  determined  values  are  shown  in  Figures  5.37  and  5.38. 

The  frequencies  are  compared  with  the  earlier  data  (Type  I 


Autocorrelation  Function 
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Figure  5.37  Strouh.il  Number 


vs. 
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of  Fiijuros  5.6  and  5.7,  rospoct  i vo ly ) . AlLhouqli  Lho  values 
of  fd/U^  obtained  from  the  v-comp<5nent  autocorrelat  ions  are 
sliijhtly  lanjer  than  those  corresponding  to  the  u-component 
crosscorre  lat  ions , they  indeed  confirm  tlie  trends  and  con- 
clusions noted  earlier. 

5 . 6 ReynoKis  Number ndependenco 

A few  runs  (of  the  Type  1,  with  symmetrically  posi- 
tioned probes)  were  taken  at  different  jet  exit  Reynolds 
nuuvber  to  determine  if  there  is  any  strong  Reynolds  number 
dependence  of  the  flopping  frequency.  Duo  to  inherent  lim- 
itations of  the  jet  setup,  only  two  different  values  were 
used. 

Measurements  of  erosseoire lat ion  functions  were 
taken  at  x/d  stations  of  *10,  (>0  and  80  for  exit  Reynolds 
numbers  of  7000  and  15100.  Fiiiure  5.30  shows  the  corre- 
sponding results  in  terms  of  tlie  dimensionless  flapping  fre- 
quency fj^d/H^.  Included  in  Pigure  5.30  are  the  results  for 
Re  = 10000  previously  reported  in  Figure  5.6. 

It  is  obvious  in  Figure  5.30,  that  tlie  frequency  of 
flapping  expressed  througli  Hi*'  parameter  f j^d  is  indepen- 
dent of  tlie  Reynolds  number,  at  least  in  the  range  investi- 
gated. This  Reynolds  numbt'r  similarity  is  a well  founded 
result.  The  flapping  bt'havior  character i'.ted  by  relatively 
low  frequencies  should  be  related  to  the  large-scale  com- 
ponents of  the  flow  field,  aiul  therefore  independent  of  the 
viscosity . 
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5 . 7 Souro(vs  of  Krror  iind  Ingcciirncios 

Since  this  investigation  is  of  an  experimental  na- 
ture, it  is  appropriate  to  discuss  the  possible  sources  of 
error  and  uncertainties  in  performing  the  reported  measure- 
ments, The  most  important  aspects  to  bo  considered  are; 
inaccuracy  in  measuring  the  mean  velocity  profiles,  error 
in  measuring  the  root-mcan-s<piare  of  the  velocity  sign<ils, 
and  uncertainty  in  ostimatinu  the  crosscorrolation  fvinc- 
tions. 

It  was  mentioned  in  Chapter  3 that  mean  velocity 
profiles  were  measured  with  a home-made  total-head  Pitot- 
tube.  This  was  so  chosen  in  order  to  avoid  uncertainties 
inherent  to  conventional  static- tubes , whose  mean  readings 
vary  strongly  with  their  size  in  turbulent  flows  (hradshaw 
[l970^ ) . On  the  other  hand,  the  selected  internal  diameter 
and  wall  thickness  v.-ero  small  enough  to  interfere  as  little 
as  possible  with  the  flow  but  still  comply  with  overrent 
standards  to  give  reliable  readings. 

Measurements  of  the  root-mean-square  values  of  tlie 
fluctuating  velocities  used  to  normalize  the  correlation 
functions,  were  done  with  a true  r.m.s.  voltmeter.  This 
type  of  instrument  theoretically  gives  very  small  normalized 
errors  if  enough  averaging  time  is  allowed  (see  for  instance 
Bendat  and  Piersol  [l97lj  ) , Large  averaging  times  v;«^re  used 
in  this  investigation  by  setting  the  time  constant  of  the 
instrument  to  high  values;  the  measured  r.m.s.  values  ox- 
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hibitca  fluctuation?;  of  no  r.\oro  than  51  of  the  e..titnatcc\  value. 

The  detornination  of  the  uncertaiutien  incurred  in 
cstiP'.atlng  a crosscorrclation  function  is  a somewhat  diffi- 
cult task.  In  the  first  place,  it  is  necessary  to  realize 
that  the  Fourier  analyzer  enployed  in  this  i nvestii-jation , 
computes  crosscorre lation  function  estimates  via  Fourier 
Transform  methods;  that  is,  a cross-.spectrum  density  func- 
tion of  the  two  sampled  variables  is  first  estimated  and 
then  the  corresponding  inverse  transfom  (the  cros?-,corrcla- 


tion  function)  is  evaluated,  as  outlined  in  Appendix  B. 

This  cross-spectrum  density  function  estimate  is 
however,  an  inconsistent  estimatca.'  (Bcndat  and  Piersol 
ri97lj).  Thus,  it  is  necessary  to  perform  a smoothing  oi'^er- 
ation  over  an  ensemble  of  estimates.  If  a large  nunrber  of 
estimates  is  used  (as  in  this  investiuatioii)  , t!\c  number  of 
statistical  degrees  of  freedom  in  the  aoproximate  chi-sguare 
distribution  of  the  samples,  greatly  increasses,  reducing  anv 
required  confidei'.ce  limit. 

On  the  other  hand  a crosscorre lation  function  esti- 
mate between  two  varial^les  x(t)  and  y(t),  defined  by: 

/t-t 


-t:^- 


X ( t) y ( t+T ) dt 


(5.29) 


for  0 £ T £ T,  or  its  digital  counterp<\rt: 


Hj,y(lm) 


1 

N-1 


M-1 

?: 

n=l 


n+r 


(5.25) 


constitutes  an  unbiased  estimate,  (tliat  is. 


the  prol;ability 
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dcnaity  fuiiction  of  the  esLiiTia,tor  is  centered  at  the  true 
value  of  the  correlation  function  R^yCi)),  since  the  ex- 
pected  value  of  (5.24)  is  given  by; 


E{R  (t)}  = E x(t)y{t+T)dt} 


•T-t 

E{x(t)y  (t+T)  } dt 

0 


= R (t) 
xy 


(5.26) 


However,  it  is  very  difficult  if  not  impossible  to 
evaluate  the  mean-square-error  of  the  correlation  function 
estimate,  as  given  by  the  variance  in  equation  (5.27): 

= 2 ^ [R,.y(T)-R^y(T)]M 

= E { R^  (t)  -R^  (t)  } 

xy  xy  ^ 


= E{- 


(T-t) 


x(u) y (u+T ) X (v) y ( v+T ) dudv  } -R^  (t) 

xy 

T-t 
r / 

E{x  (u)  y (u+t)  x(v)  y (v+T)  }dudv-R^y(T) 


(5.27) 


where  u,  v arc  integration  variables. 


Equation  (5.27)  above,  involves  a mixed  fourth  order 


moment  of  x(t)  and  y (t)  which  is  quite  difficult  to  calcu- 
late. 


It  in.iy  ho  t'.luiwn  (Kolly  ot  al.  ll‘)f)71),  D'.at  for  the 
caso  of  jointly  Ciaussian  procouson  lx{t))  and  y(t)),  tho 
abovo  variance  .ipi)roaclun} : 


(5. 20) 


wlioro  M i;.:  tlu'  band  width  of  t lu'  .in.ily/.od  niqn.iln  (aijiuimod 

to  bi'  tho  sai’.u'  for  both  r.  ianal  t?)  , and  i in  the  tiuu>  delay 

where  the  inaximuin  value  of  cor  re  la  t.  ion  occurs. 

'rhus , followimj  expi  esnion  (r5.2tt),  an  approximate 

value  for  the  normalised  iiu'an-uqua re-error  cmii  be  oldained 

for  the  pri'Si'nt  iiu-anui'i'ment  s . Foi'  the  worst  (.rases  (which 

cor ri'spondt'd  to  small  valiu's  of  x/d)  , a band-width  of 

-4 

1 Klls  and  av('ra(]inq  times  of  ')12  x 1 xlO  siu'  wi're  used. 

'I’henr f ore  at  i = 0 
m 


VarlR^^(0)l  ^ 

(0)  Ixl  0 ^x'3l2xl0 

xy 


0,02 


As  far  as  t.tu'  tinu'  delay  eha  rac  l.i'i' i s t i cs  of  the 
corn'lation  functions  is  concerned  (usihI  in  ('s  I i mat.  i inj  the 
flappiutj  f rcHiui'ncy ) , no  llu'orctical  I'rror  analysis  is 
.iViiilable  in  the  1 i t er.i  t urt' . llowi'vc'r,  the  freinu'ucy  values 
obtained  in  this  i nves  t.  i (]a  t.  ion  wvMi'  nquMt  al)  1 ('  and  consis- 
tent (particularly  with  la'spi'ct  to  f ) witliin  '10?.  for  tho 
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worst  c.isos. 


In  conclusion,  the  values  roport.etl  are  expected  to 


bo  within  tho  followimj  error  bounds: 


Variable 

F.xpected  error 

U 

± 5% 

K(0) 

t 5% 

i 10% 

^2 

-t  10% 

t 10% 

e /b 
m/ 

t 5% 

The  above  expectations  arc  best  estimates  based  on 


the  used  parameters  to  compute  lliem  as  well  as  observation 
while  performing  the  measurements. 


CllA,rTKR  6 


SL’M'iAUY  or  rrrrc'i’s  and  niroussioN 

6.1  T t\t.  t o>.Ui>.'l  { on 

Tlu>  objocLive  ot  this  i'xiH'rim<'ui  i nvost  ion 
was  to  cluu  acLeri  .!i?  tlio  l*lap[i.inv)  boliavioi-  of  a piano  taii  - 
bulont  jet.  In  this  chatitor,  tho  or.r.ont  i a I l.'alma's  of  the 
flappiinj  motion  are  svmmiari  .’.oh  and  thoii-  s i an i f i oanco  is 
qualitatively  loviowo.d  in  t lu'  context  of  t lu'  mn  aa'nt  I londs 
in  tviibulont  shear  flow  research.  It  must  be  und.'tseoK'd 
that  tlie  flappinq  motion  is  a nalviral  phenoim'isui  ••.Cv'minqly 
unrelated  to  any  Coanda- type  effects;  the  latter  .irise  from 
the  tendency  of  jets  to  attach  to  ch)si'  walls,  m sudden 
expansions . 

Tho  existence  of  some  kind  of  oruani  ;’:.it  ion  in  the 
large-scale  structure  of  the  initial  imiions  of  tuilmlent 
shear  flows,  is  ntiwadays  an  aci'cpLi'd  fact,  as.  i('view('d  in 
Chapter  3.  bepending  on  the  ii.nt  i cu  I .ir  she.ir  flow,  these 
coherent  structures  are  viewed  as  landomly  s.p.iced  vortex- 
like  patterns  moving  and  incre.ising  in  site  in  the  mo.in  flow 
direction,  and  interacting  with  I'ach  otlnm-  tl>rougl\  .\  coab's- 
cenco  process. 

Two  comments  however,  slum  Id  b('  made  it\  relation  to 
the  above  deaicription . In  the  first  place',  the  s.tatc'  of 
affairs  so  described  reforfi  to  the  lu'  u-field  ii'gioi\  of 
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turbulent  shear  flows  and  the  effect  of  this  colierence  in 
the  far-ficld  self-preserving  region,  if  any,  has  not  been 
fully  studied.  On  tlie  other  hand,  the  reported  investiga- 
tions for  the  case  of  a plane  turbulent  jet  are  rather 
scarce,  even  with  respect  to  the  itiitial  regions. 

It  may  be  argvied  in  relation  to  the  second  comment, 
that  a plane  jet  might  essentially  beliave  as  a double  mix- 
ing layer,  (with  each  one  of  its  components  the  mirror 
image  of  the  otlier)  . Although  this  conci!ption  is  geeimetri- 
cally  true  for  the  initial  region,  it  is  a rather  simplis- 
tic view  of  the  plane  jot  in  its  fully  developed  region. 

However,  the  interaction  of  the  two  layers,  both  as  they 
influence  the  potential  core  in  between,  aiul  as  t liey  come 
together  in  the  fully  merging  region,  could  account  for  some 
of  tlie  distinctive  characterist  ics  of  the  plam'  jet  even 
fa  r the r down s t ream . 

In  any  case,  it  may  be  questioned  if  hidden  in  the 
general  turbulent  structure  which  cliaracteri/.es  the  self- 
preserving region  of  a plane  jet,  there  may  be  any  kind  of 
manifi'station  of  the  series  of  events  occurring  in  the 
initial  region.  The  flapping  behavior  of  tin’  jet  may  be  one 
of  these  mani fes tations . 

6 . 2 Characteristics  of  the?  Flapping  Mot  ion 

At  this  point  if  may  be  helpful  to  summarize  the 


essential  characteristics  of  such  behavior,  as  concluded 
based  in  this  investigation. 


Beginning  at  cross-sections  at  about  20  times  the 
slot  width,  crosscorrelations  taken  between  u-component  vel- 
ocity fluctuations  at  points  on  each  side  of  the  jet,  pre- 
sented a distinctive  negative  value  at  zero  time  delay.  As 
the  time  delay  is  increased,  the  correlation  functions 
assume  a quasi-periodic  shape  with  an  approximately  defin- 
able frequency.  This  result  indicates  that  both  velocity 
signals  behave  (in  the  average)  as  if  they  included  a low 
frequency  wave  component  in  antiphase  with  respect  to  each 
other.  This  could  be  the  consequence  of  the  whole  flow 
field  moving  (again  in  the  average)  in  a pseudo-oscillatory 
way. 

The  various  half-periods  between  alternative  posi- 
tive and  negative  peaks  detected  in  the  correlation  func- 
tions, are  generally  not  equal.  No  explanation  is  given  in 
this  work  for  this.  It  may  be  speculated  however,  that  in- 
asmuch as  the  flapping  motion  is  characterized  by  a rela- 
tively low  frequency,  it  may  be  coupled  with  the  low- 
frequency  components  of  the  turbulent  velocity  field.  Any 
degree  of  correlation  between  such  components  of  the  two 
measurement  points  in  the  flow  field  might  in  turn  be 
affected  by  the  randomness  of  the  high-frequency  components. 
As  the  time  delay  between  the  signals  increases  from  the  one 
corresponding  to  the  highest  degree  of  correlation,  the 
two  signals  should  become  less  dependent  on  each  other.  The 
quasi-periodic  nature  of  the  correlation  function 
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(attributed  to  the  flapping  motion)  should  then  become  less 
defined.  Recall  that  only  for  perfectly  periodic  signals 
would  the  correlation  between  them  bo  also  perfectly  peri- 
odic. 

Support  to  the  above  may  bo  derived  by  examining  the 
data  plotted  in  Figure  6.1.  It  was  obtained  by  first  fil- 
tering both  signals  before  processing  them.  Low-pass  fil- 
ters, with  cut-off  frequency  of  25  Hz,  were  used.  The 
periodic  yet  decaying  nature  of  the  correlation  function  is 
noted. 

6 . 3 Possibility  of  Puffing 

During  the  earlier  stages  of  this  investigation,  the 
existence  of  "puffing"  was  raised  as  the  possible  cause  for 
the  difference  between  the  various  half-periods  in  the  cor- 
relation functions  (and  possibly  the  decaying  correlation) . 
This  puffing  effect  might  be  originated  by  instabilities  in 
the  blower  feeding  the  nozzle.  Correlation  functions  were 
ta)cen  with  various  probe  arrangements  in  the  potential  core 
close  to  the  nozzle  exit  (see  example  in  Figure  6.2).  Any 
puffing  would  have  to  be  accompanied  by  a periodic  correla- 
tion, positive  near  x = 0.  These  tests  revealed  no  abnor- 
malities which  could  indicate  the  existence  of  this  effect. 

6 . 4 Frequency  of  Flapping  Motion 

The  frequency  of  flapping,  although  almost  constant 
at  each  cross-section,  decreases  in  the  downstream  direction. 
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In  the  far  region  it  seems  to  scale  ijuite  well  by  the  in- 
verse of  a characteristic  development  time  of  tlie  mean 
flow,  given  by  the  ratio  of  the  half-widtli  to  the  moan  vel- 
ocity at  the  centerline  of  the  jet. 

The  amplitude  of  the  correlation  function  at  zero 
time  delay  R(0) , depends  on  the  lateral  coordinate.  For 
small  separations  it  definitely  should,  since  it  corre- 
sponds to  a space-correlation  for  varying  separation  dis- 
tance between  probes  (equation  5.8,  Figure  5.2*1).  These  re- 
sults when  compared  with  the  ones  reported  by  t'utmark  and 
Wygnanski  (1976)  showed  some  disagreement  (see  Figure  5.24). 
Fven  though  the  present  work  shows  the  same  trends  as 
Gutnark  ami  Wygnanski  [1976]  for  space-correlat ions  corre- 
sponding to  y/b  £ 0.5,  tliere  is  noticeable  disagreement  with 
the  cases  at  y/b  2.  1*0.  Particularly,  the  curv'e  correspond- 
ing tc  y/b  = 1.5  reported  by  the  above  mentioned  aut’aors, 
shows  a negative  dip  at  ajjpro.x imate ly  Ay/b  = 2.3  with  a 
trend  to  become  more  negative  as  the  abscissa  increases. 

This  is  an  unlikely  result  since  it  should  be  realized  that 
having  the  probes  separated  by  Ay/b  = 2.3  witli  middle-point 
at  y/b  =1-5  corresponds  to  effectively  setting  one  probe 
at  y/b  = 0.35  and  the  second  one  at  y/b  = 2.65.  The  latter 
probe  is  in  the  entrainment  region  totally  ovitside  of  the 
turbulent  core,  and  any  correlation  between  t)ie  x-component 
of  the  velocities  at  tliese  points  sliould  be  close  to  zero. 
That  is  not  the  case  with  Gutmark  and  Wygnanski  (1976]  but 
is  indeed  so  in  the  present  work. 
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6.5  Coi  voction  of  Flat>pinq  Motion 

Correlation  measurements  with  the  probes  at  differ- 
ent x/d  stations  and  yet  on  opposite  sides  of  the  jet  showed 
that  the  flapping  motion  seems  to  travel  downstream.  This 
flapping  convection  velocity  roughly  corresponds  to  the  tur- 
bulent convective  velocity  previously  measured  and  reported 
by  other  investigators  (Ott  (1972),  Young  [1973]). 

6. 6 Relationship  to  Wave-Guide  Theory 

In  trying  to  put  the  described  results  in  a proper 
theoretical  framework,  one  is  tempted  to  consider  the  cur- 
rent efforts  in  the  stability  theory  of  turbulent  shear 
flows.  It  should  be  stated  however  that  no  mathematical 
modeling  is  intended  here,  but  instead  a sununary  of  the  main 
features  of  such  theory  in  the  light  of  the  experiments 
reported  in  this  work.  A review  of  spatial  stability  theory 
is  given  in  Appendix  C. 

A realistic  wave-guide  model  of  turbulent  shear 
flows  should  account  for  a spatially  growing  wave  of  given 
frequency,  attaining  its  maximum  amplitude  through  some  kind 
of  saturation  effect  and  decaying  in  the  moan  flow  direction. 
In  fact,  as  Liu  [1974]  has  suggested,  it  is  possible  to  con- 
struct such  a wave-gviide  representation  by  assuming  the  tur- 
bulent flow  to  be  composed  of  three  parts:  the  mean  flow, 
the  periodic  wave  component  and  the  small-scale  turbulence. 
The  organized  wave  component  may  be  assumed  as  the  product 
of  an  amplitude  function  to  be  determined  from  the  kinetic 


I 
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energy  integral  equation  and  a shape  function  obtained  from 
the  linear  eigenvalue  problem  corresponding  to  the  local 
mean  flow.  The  various  flow  quantities  (assumed  to  include 
such  wave  component)  are  introduced  into  the  basic  equation 
of  motion  and  time-  and  phase-averages  are  performed.  An 
integral  equation  can  in  turn  be  obtained  for  the  kinetic 
energy  of  the  instability  wave  and  numerical  solutions  can 
be  sought  after  considering  proper  boundary  conditions. 

Throe  essential  characteristics  for  the  organized 
wave  are  obtained  through  this  type  of  analysis  (Merkino  and 
Liu  [1975]):  (a)  the  amplitude  of  the  wave  increases,  at- 

tains its  maximum  value  and  tlien  decays  in  the  downstream 
direction,  mainly  through  kinetic  energy  exchange  mechan- 
isms; (b)  the  lower  the  frequency  of  tlie  wave  component,  the 
further  downstream  its  maximum  amplitude  is  attained; 

(c)  the  sinuous  modes  (that  is,  the  organized  waves  oscil- 
lating in  antiphase,  like  the  u-componont  velocities  in  the 
flapping  behavior)  have  larger  amplitude  rates  and  show  more 
delayed  saturation  values  in  the  downstream  direction,  thv\n 
the  varicose  (or  symmetrical)  modes. 

It  is  seen  therefore,  that  a wave-guide  representa- 
tion could  bo  used  in  modeling  some  of  the  experimental  ob- 
servations reported  in  this  work.  For  instance,  the  sinuous 
modes  (the  most  amplified  ones)  predicted  by  stability 
theory  might  account  for  the  pseudo-oscillatory  u-components 
of  the  velocity  (in  antipliase)  cliaracter izing  the  flapping 


motion.  On  the  other  hand,  the  theoretical  result  that  the 
lower  frequency  wave  components  have  the  tendency  to  reach 
their  maximum  growth  further  downstream,  might  be  related 
to  the  relatively  low  and  decreasing  frequencies  associated 
to  the  flapping  behavior. 

In  concluding  this  chapter,  a perspective  of  the 
plane  turbulent  jet,  might  be  as  follows.  In  the  initial 
region,  the  mixing  layers,  one  at  each  side  of  the  potential 
core,  seem  to  be  characterized  by  the  presence  of  large- 
scale  coherent  structures.  As  the  two  layers  merge  togeth- 
er, the  coherent  structures  interact  through  a yet  unknown 
process.  As  a result,  the  possible  existence  of  organized 
waves  as  predicted  by  stability  theory,  might  bo  manifested 
in  the  far  downstream  self-proserving  region  in  terms  of 
distinctive  phenomena,  among  them  the  flapping  motion. 


CHAPTER  7 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  WORK 

The  following  conclusions  cun  bo  ostublished,  bused 
on  the  experimon tu 1 results  obtained  in  this  investigation: 

1)  The  flapping  motion  of  a turbulent  plane  jet  is 
a distinctive  and  measurable  natural  phenomenon.  This  was 
confirmed  by  spot-checks  at  two  other  jet  setups  and  by  a 
previous  preliminary  investigation  (Goldschmidt  and  Brudshav/ 
119731) . 

2)  The  frequency  of  flapping  decreases  in  the  longi- 
tudinal direction  and  remains  unchanged  in  the  lateral  di- 
rection . 

3)  Approximate  self-preservation  is  obtained  for  the 
flapping  frequency  (for  x/d  > 30),  if  scaled  with  the  cen- 
terline mean  velocity  and  the  half-width  of  the  jet,  giving 

fb/U^  = 0.11. 
m 

4)  An  estimate  of  the  amplitude  of  flapping  gave 
values  in  the  order  of  20%  of  the  jet  half-width. 

5)  The  flapping  behavior  travels  in  the  downstream 
direction  with  a velocity  about  25%  smaller  than  the  turbu- 
lent structure  convective  velocity. 

6)  The  frequency  of  flapping  is  independent  of  the 
Reynolds  number  in  the  range  7900  ^ Rc  £ 15100. 
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Further  work  is  recommended  in  the  follow! ncj  areas: 

1)  More  accurate  measurements  are  needed  for  x/d  < 
20.  The  use  of  a nozzle  with  a larger  slot  would  permit 
bettor  accuracy  in  positioning  the  probes  at  those  stations. 

2)  Simultaneous  data  acquisition  in  more  than  two 
points  of  the  flow  field,  can  be  accomplished  by  using 
magnetic  tape  recording  instrumentation.  This,  together 
with  a large  computer,  would  increase  the  versatility  of 
data  reduction  and  analysis. 

3)  A detailed  and  accurate  experimental  study  is 
suggested  for  the  initial  zone  of  the  turbulent  plane  jet. 
Particularly  in  the  region  where  the  two  mixing  layers 
merge  together.  The  aim  of  this  investigation  would  bo  to 
get  a better  insight  on  how  they  interact  with  each  other. 

4)  Preliminary  and  limited  observations  with  an  ex- 
ternal acoustic  input  (frequencies  less  than  600  liz)  showed 
no  changes  in  the  flapping  motion  of  the  jet.  More  experi- 
ments are  needed  in  this  respect. 

5)  Mathematical  modeling  of  turbulent  shear  flows 
through  stability  theory  (although  in  a developing  stage)  is 
a promising  one.  A continuous  effort  in  this  area  and  its 
possible  relationship  with  the  flapping  behavior  of  a tur- 
bulent plane  jet,  is  strongly  recommended. 
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Appendix  A 

Self-Preservation  Character ist ics 
of  the  Plane  Turbulent  Jet 
Mean  velocity  traverses  were  performed  at  several 
x/d  stations,  as  described  in  Chapter  3.  The  resulting  pro- 
files are  presented  in  Figure  A.l,  where  the  local  mean 
velocity  and  the  lateral  y-coordinate  have  been  nondimen- 
sionalized  by  the  centerline  mean  velocity  U and  the  jet 
half-width  b respectively,  corresponding  to  each  x/d  sta- 
tion. Also  presented  in  Figure  A.l  are  Reichardt's  (1951] 
experimental  results.  Self-preservation  of  the  turbulent 
jet  is  then  attained  for  the  range  of  x/d  reported  in  this 
i nvestigation. 

The  similarity  and  self-preservation  principles, 
already  well-documented  both  theoretically  and  experimental- 
ly in  the  existing  literature  (see  for  instance,  Abramovich 
[19G3],  Schlichting  [1968]),  allow  one  to  write  the  follow- 
ing e.xpressions : 


b 
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"1  - ^'1^ 


(A.  1) 


U ^ 
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In  these  relationships,  Kj^,  K2/  C2 
constants  the  first  two  measuring  the  widening 
flow  and  the  last  two  indicating  the  geometric 
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virtual  origins  of  the  jet,  respectively.  There  is  no  rea- 
son to  expect  coincidence  of  those  origins  as  reported  by 
Flora  and  Goldschmidt  [1969].  Being  their  determination 
the  result  of  graphical  extrapolation  and  therefore  not 
sufficiently  accurate,  their  values  are  different  from  one 
experimental  setup  to  another  as  noted  in  Table  A.l,  where 
those  parameters  obtained  by  various  authors  are  compared. 

The  results  obtained  in  the  present  investigation 
corresponding  to  expressions  (A.l)  and  (A. 2)  respectively, 
are  as  follows: 


I = 0.083  (j  + 6.62) 


{A.  3) 


U 

(gSi)  = 0.24  (J  -4.53) 


(A.  4) 


They  are  also  presented  graphically  in  Figures  A. 2 and  A. 3, 
respectively. 


Table  A.l  Spreading  Parameters  of  Plane  Jets. 


Ci.T.ensionless  Centerline  Mean  Velocity  (U 


i X B 

App  1 ic.t  t i of  t ho  HP  Fouf  icT  Aiui  1 y /or 

Tlu’  Ik'wli'f  t,-Packoi  d 54  52A- I 4 H Fourior  Analy/.or 
dystom  is  able  to  com[)iito  the  crosscorrelation  fnnction  be- 
tween any  two  given  analoc  voltage  signals.  noth  signals 
representing  any  two  physical  cinantities,  are  first  low- 
pass  filtered  (accordii\g  to  the  sampling  r.itt'  useil  in  con- 
verting the  signals  to  digital  form)  in  ortler  to  avoid 
aliasing.  Once  tligiti/.ed  and  stored  in  the  memory,  the 
crosscorrelation  between  tlie  signals  is  competed  via  the 
Fast  Fourier  Transform  algorithm,  as  t lie  flow  diagram  in 
F i g u re  It  .1  i nd  i ca  te s . 

Th('  vinaly/er  can  do  this  operation  in  a repetitive 
way  for  any  number  of  samples  and  com[UiLo  an  .iveianje  of  t lie 
if  properly  programmed.  The  block  diagram  and  listing  of 
a program  used  to  estimate  an  averagi'  of  crosscorre l.it  ion 
functions  from  250  samples,  are  [irescsited  in  Figure  It . 2 aiul 
Table  It.  1 , respectively. 

Tn  order  to  estimate  the  aut ocor rel a t i on  function 
of  the  v-coniponent  of  the  velocity  at  the  centerline,  a dif 
ferent  and  somehow  more  complicated  proqram  was  used.  In 
fact,  it  was  nece.ss.iry  to  compute  actually  four  different 
correlation  functions  for  this  case,  as  egiuition  (It.l) 


statc.s : 
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Figure  B.l 

Block  Diagram  of  Crosscorrola tion  Computation 


I 
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Table  B.l  Program  Ui;od  by  the  Fourier  Analyzer  in  Computing 
Crosscorrelation  Averages . 


LABEL  0 ENTER 
CLEAR  2 ENTER 
LABEL  1 ENTER 

ANALOG  0 SPACE  1 SPACE  2 ENTER 

CLEAR  1 SPACE  0 SPACE  128  ENTER 

CLEAR  1 SPACE  3G4  SPACE  512  ENTER 

CORR  1 ION  TER 

+ 2 ENTER 

STORE  2 ENTER 

COUNT  I SPACE  250  ENTER 

f 0 SPACE  250  ENTER 


END  ENTER 
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R^(t)  = V(t)v(t+T)  = 


4S' 


(t)  (t+T  ) + e2(t)e2(t+T) 


- (t)  02  (t+T  ) - eJUTe^Tt+rT 


(B.l) 

where  and  are  the  voltage  signals  from  the  two 

sensing  elements  of  the  x-wire,  respectively,  and  is  the 
v-component  sensitivity  of  the  wires.  (S^  is  approximately 
the  same  for  both  wires  if  they  are  properly  matched) . 
Equation  (B.l)  can  be  obtained  easily  with  some  algebraic 
manipulations,  from  the  standard  relationships  defining  the 
voltages  of  an  x-wire: 


Of (t)  = u (t)  + V (t) 


e2(t)  = u(t)  - v(t) 


(B.2) 


A listing  of  the  program  used  for  the  autocorrela- 
tion computation  is  given  in  Table  B.2. 

In  estimating  any  function  which  involves  time  shift- 
ing operations  (like  correlations  or  convolutions) , the 
analyzer  computes  integrations  between  the  limits  0 and  T 
(that  is,  the  size  of  the  sampled  data).  The  analyzer 
further  assumes  that  the  sampled  data  is  periodic  with 


It 
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Table  D.2  Program  Ur.od  by  the  Fourier  Analyr-or  in  Computing 
Autocorrelation  Averages • 


LABKL  0 KUTER 
CLEAR  3 ENTER 
LABEL  1 ENTER 

ANALOG  4 SPACE  5 SPACE  3 ENTER 

LOAD  4 ENTER 

STORE  1 ENTER 

SUBJUMP  2 ENTER 

STORE  2 ENTER 

LOAD  5 ENTER 

CORR  1 ENTER 

INTERCUNG  2 ENTER 

- 2 ENTER 
STORE  2 ENTER 
LOAD  5 ENTER 
STORf:  1 ENTER 
SUBJUMP  2 ENTER 
+ 2 ENTER 
STORE  2 ENTER 
LOAD  4 ENTER 
CORR  1 ENTER 
INTERCUNG  2 ENTER 

- 2 ENTER 
+ 3 ENTER 


STORE  3 ENTER 


Table  B. 


IIP  M - 
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(continued) 


COUNT  1 SPACE  n ENTER 
v 0 SPP.CE  n ENTER 

— f 

SUBROUTINE 

CLEAR  1 SPACE  0 SPACE  128  ENTER 
CLEAR  1 SPACE  384  SPACE  512  ENTER 
CORR  1 ENTER 


SUBRTRN  ENTER 
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period  T.  The  estimated  function  therefore  results  in  error 
(the  so  called  wrap-around  error).  In  order  to  effectively 
avoid  this  (at  least  partially) , one  of  the  data  records  is 
cleared  to  zero  from  0 to  T/4  and  from  3/4T  to  T,  before 
the  data  is  processed.  In  this  way,  the  resulting  cross- 
correlation function  is  valid  from  0 to  T/4  and  from  3/4T 
to  T (corresponding  to  the  interval  -T/4  £ t ^ T/4).  The 
remaining  parts  of  the  estimated  function  are  still  in  error 
(IIP  Fourier  Analyzer  Training  Manual)  . 
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Soino  Aspect.:;  of  St.ilMHty  'L'liu o r y 

Considerable  interest  has  developed  in  modeling  the 
large-scale  structure  of  turbulent  rliear  flows  through  sta- 
bility analyses.  This  [parallel  aptiroiich  has  Loon  mainly  mo- 
tivated by  the  jet-noise  geiieraticn  [Jichlem. 

The  objective  of  any  theoj'y  of  hyilrijilyiiaini  c fitiibi  l- 
ity  is  to  predict  the  future  behavior  of  a viisturbance , once 
it  is  introduced  into  the  flow  being  analyr-ed.  In  tlae 
classical  linear  theory  of  parallel  flows,  small  perturba- 
tions are  considered  in  the  equations  of  motion  and  these 
in  turn  arc  linearized. 

Assuming  that  any  disturbance  is  composed  of  travel- 
ling waves  witli  an  exponential  time  factor,  the  equations 
of  motion  for  the  case  of  a two-dimensional  flow  with  two- 
dimensional  disturbances,  reduce  to  a fourth-order  differen- 
tial equation  (the  so-called  Orr-Son'imerfc  Id  ecpiation)  . This 
equation  togetlier  with  appropi'iate  bovindary  conditions,  con- 
stitutes an  eigenvalvie  problem,  the  solution  of  which,  de- 
termines the  oscillatory  modes  v;iLh  th.eir  clia  racte  r i s t ic 
wave-nun\ber,  propagation  speed  and  amp  li  f ication  factor.  If 
solutions  are  found  that  decay  in  time,  th>;  flow  is  consid- 
ered as  stable;  if  they  increase  with  time,  tlie  flow  is 
said  to  be  unstable  (Sch licliting  [I'K^S]). 


The  theory  has  been  applied  witli  considerable  suc- 
cess to  some  simple  laminar  flows  in  the  verge  of  Ijccoming 


tuxbulcat,  .uid  the  swnilAtity  between  the  turbulent  coherent 
structures  and  the  events  leadin^j  to  transition  in  laminar 
flows,  liave  motivated  some  investigators  to  study  the  former 
through  stability  theory. 

hydrodynamic  stability  tlieory,  as  applied  to 
laminar  and  transitional  regimes  of  mixing  layers  and  jots, 
constitutes  a subject  by  itself.  /'juong  the  mutable  works 
in  this  area  are  the  ones  by  MichalVve  ^f9  70j  (for  \i\stai\ce)  , 
Mattingly  and  Criminalo  [l97l],  Batclielor  and  Gill  [1962], 
and  Caster  [1965]  covering  the  theoretical  aspects,  and  the 
ones  by  Sato  [i960] , Phillips  [1966] , Froymuth  [1966]  and 
Miskad  [1972]  dealing  with  physical  experiments. 

T}ie  work  by  Landahl  |1967]  is  an  example  of  this 
postulated  relationship  between  the  turbulent  structure  and 
the  events  related  to  transition.  Tn  attempting  to  relate 
the  statistical  properties  of  tlie  wall  pressure  in  a turbu- 
lent boundary  layer  to  the  characteristics  of  the  mean  flow, 
Landahl  [l96?]  formulated  a wavc-yuidc  miOdel  (i.e.,  the 
admission  of  wave  propagation  modes)  for  the  turbulent  vol- 
ocity  f luctua ti ons . Aii  aigjroximatc  solution  v/.is  obtained 
from  a non-homogeneous  Orr-Sommor fe Id  cquatiori,  where  fhe 
non-linear  turbulent  stresses  were  considovod  as  the  forcing 
elemciits.  Consiilorablo  avjrcement  with  tlio  cxpcrimoiita  1 re- 
sults reported  by  Corcos  [196 -ij  was  obtained  for  tlu’  stream- 
wisG  decay  of  t)ie  fluctuations.  However,  tin}  aejrccmcnt  was 
not  found  for  the  measured  and  predicted  convective  veloci- 


One  method  of  getting  an  insiglit  into  shear  flow 
structure  and  experimentally  testing  v/ave-guide  models  is  to 
introduce  a periodic  disturbance  of  known  frequency  and  am- 
plitude into  the  flow  and  to  determine  the  response  of  the 
flow  to  it.  Such  a procedure  was  followed  by  Crov/  and 
Champagne  [l97l]  for  a circular  jet.  They  introduced  a per- 
iodic surging  upstream  of  the  plenum  chamber  and  measured 
the  response  along  the  centerline  in  terms  of  several  axial 
component  quantities  like  mean  velocity,  root-mean-square 
fluctuations,  etc.  After  trying  various  frequencies  and 
eunplitudes,  they  obtained  as  a preferred  mode,  the  one  cor- 
responding to  the  average  frequency  of  vortex-puffs  formation 
in  the  transitional  region  of  the  jet,  previously  measured 
by  themselves  through  visualization  experiments  (it  corres- 
ponded to  a Strouhal  nurriber  fd/U^ofO.3).  According  to 
their  results,  the  preferred  wave  attained  its  maximum  am- 
plitude through  the  combined  effects  of  linear  amplification 
and  non-linear  saturation  at  a certain  axial  position  and 
then gradua lly  decayed  downstream. 

Following  the  same  pattern,  Chan  [1974]  performed 
pressure  measurements  in  a disturbed  circular  jet.  Through 
an  appropriate  triggering  scheme  using  the  narrov/-band 
filtered  signal  from  a traversing  microphone,  ho  obtained 
oscilloscope  traces  representing  the  spatial  variation  of 
the  pressure  signal  at  a fixed  time.  The  results  obtained 
along  the  centerline  agreed  with  the  ones  by  Crow  and 


ChampacjriG  i:,nd  Ko  aiul  DAYicis  [1971].  Furthermore,  his  mee- 
suremc.uts  for  the  mixing  layer  and  the  near  field  seemed  to 
confirm  a stability  model  previously  dovelo[jcd  by  Michalke 
[1972],  as  referred  to  by  Chan  [1974]. 

In  getting  appropriate  data  for  the  development  of 
the  wave-guide  approach  of  shear  flow  turbulence,  Hussain 
and  Reynolds  [ 1970,  1972  ] perfonned  experiments  on  a 
channel  flow  with  artificially  introduced  wave  disturbances. 
Using  a phase  averaging  tech.nique  with  selective  sampling 
through  the  cyclic  position  of  the  wave-maker,  they  surveyed 
downstream  of  the  initial  disturbance  and  extracted  the  or- 
ganized wave  motion  from  a hot-wire  signal.  Realizing  that 
the  experimental  results  suggested  that  more  than  one  wave 
mode  was  excited  by  the  initial  disturbance,  they  performed 
an  approximate  single-mode  linear  analysis  of  the  propagation 
characteristics.  The  dynamical  equations  were  further  de- 
rived for  the  small  amplitude  wave  disturbances  and  closure 
schemes  were  proposed. 

More  recently,  Hussain  and  Zaman  [1975]  reported  pre- 
liminary results  in  extending  the  previous  experiments  by 
Crow  and  Champagne  [l97l]  to  a wider  range  of  excitation 
Strouhal  numbers. 

The  instability  waves  developing  in  compressible 
flows  v/ere  mathem.atically  studied  by  Liu  [1974]  for  a mixing 
layer  and  by  Merkine  and  Liu  [1975]  for  a plane  jet.  They 
considered  the  kinetic  energy  aspects  of  the  organized  wave, 
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once  it  Kns  i>cci\  sopc\iriXtcJ  iiron"\  the.  ir^cnn  nncl  iinc^scnlc  corn." 
ponents  of  tlio  flow.  A remark  wan  made  on  tlie  role  of  the 
synvuotric  or  variconc  and  asynunetric  or  yinuous  modes. 

The  v/avo-guide  representation  of  turbulent  shear 
flows  has  some  inherent  difficulties.  For  instance,  the 
case  of  interest  in  a physical  problem  is  included  in  spa- 
tial stability  theory.  However,  this  foimulation  predicts 
an  exponential  amplification  of  the  wave  in  the  downstream 
direction  and  the  original  small  amplitude  assumption  be- 
comes inconsistent,  unless  some  hind  of  limitation  is  im- 
posed on  that  amplification. 

In  any  case,  a realistic  wave-guide  model  should 
account  for  a spatially  growing  wave,  v/nich  attains  a maxi- 
mum growth  through  some  kind  of  saturation  effect  after  non- 
linear interactions.  Thereafter  it  might  begin  to  decay, 
as  the  experiments  of  Crow  and  Champagne  [lOVlJ  suggested 
for  the  circular  jet. 

One  way  of  studying  turbulent  shear  flows  through  a 
wave-guide  representation,  is  by  separating  the  flow  field 
into  three  components:  the  mean  flov/,  the  instability  wave 
and  the  small-scale  turbulent  fluctuations.  Once  the  fluctu- 
ating quantities  liave  been  decomposed  in  the  above  three 
terms,  they  are  then  introduced  into  the  basic  equations. 

An  ordinary  tim.e  average  is  taken  first,  followed  by  a phase 
averaging.  This  leads  to  dynamical  equations  for  each  one 
of  the  components.  The  phase  averaging  is  defined  as  the 


cworAye  over  a l^r'jo  nunUjcr  o£  porioilu  of  the  \y\\\Q  at  a 
ijivcn  position  (Liu  [1974])  and  is  physically  realizable 
for  instance,  by  considerintj  an  cnseiiilile  of  j)ointG  having 
the  same  phase  with  respect  to  a reference  oscillator,  wlicn 
a disturbance  is  artificially  introduced  (Hussain  and 
Reynolds  [l9  70]). 

Reynolds  and  Ilusr.ain  [19  72]  pointed  out  a serious 
closure  problem  ii\  the  equatioi\s  for  the  orvia:i i zed  wave  as 
the  result  of  ttie  outlined  averaging  procedures.  They 
further  considered  the  necessity  to  take  info  account  the 
interaction  of  the  organized  waws  aiul  the  small  scale 
f luctuations , perhaps  through  the  use  of  an  e<idy  viscosity 
model,  as  Landahl  [1967]  ha^-l  prcviovisly  sjioculated. 

Following  the  same  aj^proach  of  scipara ti nij  the  flow 
field  into  three  components,  Liu  [1974]  formulated  a mode  1 
for  a compressible  mixing  l.iyer  and  later  c?xtendcd  it  to  a 
two-dimensional  jot  (?lerkine  and  Liu  [1979]  ) . in  wliat  fol- 
lows, a sunmiary  of  such  procedure  is  given. 

Any  flow  quanta' ty  wliich  is  a function  of  time  .ind 
spatial  coordinates  is  split  into  three  components: 

q = q + q ' + q " ( C . 1 ) 

whore  q denotes  the  time-ave raged  compoiient,  ii  ' the  ii,  - 
stantaneous  periodic  compoiient  and  q"  the  siita  1 1 -sea  1 e turbu- 
lent componor.t.  Tlio  organized  w.ave  pa.rt  q'  is  given  by: 


q’  = < q > - q 


(C.2) 


where  < > indicates  phase  averaging. 

Introduc  ng  these  quantities  into  the  Navier-Stokes 
equations,  performing  the  described  averaging  procedures 
and  making  use  of  boundary  layer  type  of  approximation,  a 
two-dimensional  von  Kdrmdn  integral  equation  can  be  ob- 
tained for  the  kinetic  energy  of  the  instability  wave,  which 
for  an  incompressible  flow,  roads: 


1 A 

2 dx 


u (u ' ^+v ' ^ ) dy  = - 


u'v'  A dy 


(<u"^>-u"^) 

J X 


(<u"v">-u"v")  (A'4.  A'  ) 
9y  3X 


2 2 V ' 

-(<iy"  >-v"  ) ~ dy 


(C.  3) 


The  left-hand  side  of  equation  (C.3)  represents  the 
net  rate  of  increase  of  the  kinetic  energy  of  the  wave  and 
the  terms  on  the  right-hand  side  arc  the  mechanisms  govern- 
ing this  change:  energy  conversion  from  the  mean  flow  (or 
production)  and  kinetic  energy  exchange  between  the  wave  and 
the  fine-scale  turbulence  (molecular  viscosity  effects  have 
been  neglected) . 
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If  it  is  further  assumed  a linear  relationship  be- 
tween the  wave  induced  stresses  and  the  instability  wave 
rate  of  strain  of  the  type: 

-(<u"^>-u"^)  = 2e|^’  (C.4) 

then  equation  (C.3)  reduces  to: 
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(C.5) 


where  e is  an  eddy  viscosity. 

Once  the  basic  equation  has  been  established,  clo- 
sure is  attempted  (and  here  is  the  essence  of  Liu's  model) 
by  assuming  the  wave  as  the  product  of  an  amplitude  function 
(to  be  determined  from  the  kinetic  energy  equation) , and  a 
shape  function  (to  be  determined  from  the  linear  eigenvalue 
problem  corresponding  to  the  local  moan  flow) . 

In  effect,  any  fluctuating  quantity  q'  pctaining  to 
the  organized  wave  is  assumed  to  be  given  by: 


q'(C,n)  = A(^)Q(n;fJcxp(-i3t)  + c.c. 


(C.6) 
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where  Q is  the  ciycnf unction  obtained  from  the  local 
spatial  stability  theory,  A is  the  modifying  amplitude 
function,  and  F,  = x/d,  n - y/b,  0 = local  dimensionless 
frequency,  t = dimensionless  time,  and  c.c.  stands  for 
complex  conjugate. 

The  eigenvalue  problem  to  determine  the  shape  func- 
tion Q is  solved  with  appropriate  boundary  conditions- 
Among  them,  the  possibility  to  have  varicose  or  synmetric 
and  sinuous  or  asymmetric  (the  one  connected  to  the  flapping 
behavior)  modes  are  considered. 

As  discussed  by  Merkine  and  Liu  [l975]  the  amplifi- 
cation and  decay  of  the  organized  wave  comes  primarily  from 
the  amplitude  function  A rather  than  from  the  amplification 
rates  of  the  stability  theory.  The  eigenfunctions  of  the 
type  Q and  their  associated  amplification  rates,  play  only 
a subsidiary  role. 

The  amplitude  function  may  bo  determined  by  substi- 
tuting the  equations  of  the  type  (C.6)  into  the  energy  re- 
lation (C.5).  An  expression  for  the  amplitude  function  is 
obtained  involving  the  various  energy  integrals;  those  can 
be  evaluated  numerically  after  assuming  an  initial  condition 
through  the  energy  do.  ■ ity  of  the  initial  disturbance. 

Although  the  computational  results  obtained  by 
Kerkine  and  Liu  are  not  directly  a.ppliciible  to  the  present 
investigation  (as  they  refer  to  a high-speed  compressible 
jet  flow),  some  of  their  findings  may  be  cjual i tati ve ly 
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described. 

The  Ycxrious  enorejy  exchr^ngo  mochani ;imr.  in  etiuation 
(C.5)  determine  the  development  of  the  amplitude  of  the  or- 
yani/cd  wave.  In  the  first  stayes,  the  amplitude  increases 
rapidly  as  the  result  of  laryor  production  ;;s  compared  with 
dissipation.  As  the  development  continu.^s,  the  dissipation 
term  dominates  over  the  production  one,  causiny  tlie  ampli- 
tude of  the  wave  to  attain  its  maximum  value  an.d  t!\en  to 
decay . 

The  above  described  process  is  retarded  for  the  lew 
frequency  components;  in  other  words,  the  lower  the  fre- 
quency of  the  wave,  the  further  downstream  its  maximum  am- 
plitude is  attained. 

Finally,  the  role  of  the  varicose  and  simuous  modes 
(althouyh  not  much  differentiated  at  tlie  early  stayes)  is 
such  that  the  sinuous  waves  fiave  laryer  a;nf>]  i f i ca  tion  rates 
(as  predicted  by  local  stability  theory  and  suyyested  by 
earlier  workers  (Rockwell  and  Niccolls  [l972l  , Mattinnly  and 
Criminalc  flOTl]))  and  show  more  delayed  saturation  values 
in  the  downstream  direction. 

The  formulation  of  a wave-yuide  theory  for  turbulent 
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shear  flows  is  a very  promising  one,  but  is  still  in  a de- 
veloping st.aye  with  many  difficulties  to  he  overcome. 
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